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Abstract 


This  Final  Technical  Report  describes  completed  research  accomplishments  and 
ongoing  activities  that  are  focused  on  the  evolution  of  boundary  driven  acoustic  dis¬ 
turbances  in  a  low  Mach  number  shear  flow  like  that  found  in  the  chamber  of  a  solid 
rocket  engine.  The  completed  work  (manuscripts  in  Appendices  A  and  B)  focuses  on 
the  relatively  complex  wave  systems  that  appear  in  a  two-dimensional  planar  shear 
flow  following  the  refraction  of  very  simple,  initially  planar  axitil  disturbances.  Work 
m  progress  emphasizes;  (1)  the  characteristics  of  acoustic  disturb«mces  driven  by  side- 
wall  mass  addition  in  semi-confined  channels  and  tubes,  (2)  the  role  of  “strongly  in¬ 
jected”  Stokes  boundary  layers  in  providing  a  transition  from  the  acoustic  flow  to  the 
no-slip  condition  on  the  wall,  and  (3)  mathematical  methods  required  to  deal  with 
nonlinear  processes  within  an  acoustically  disturbed  flow.  The  review  of  our  work 
emphasizes  the  importance  of  studying  the  evolution  of  boundary  driven  acoustic  dis¬ 
turbances,  primarily  to  gaun  an  understanding  of  how  small  burning  rate  trainsients 
(modelled  by  unsteady  wall  injection)  lead  to  large  engine  chamber  responses  observed 
in  unstable  solid  rockets. 


Chapter  1 
Introduction 


Burning  rate  transients  of  a  combusting  solid  propellant  produce  gaseous  products 
that  are  “injected”  into  the  rocket  engine  chamber  as  ctn  unsteady  flow.  This  type 
of  transient  mass  injection  is  known  (Kassoy  1979)  to  be  a  source  of  mechanical  aind 
thermodynamic  disturbances  that  can  be  described  in  terms  of  acoustic  phenomena. 
Our  research  program  is  focused  on  relating  the  evolution  of  acoustic  phenomena  in 
semi-confined  low  Mach  number  shear  flows  to  specific  types  of  boimdary  driving, 
particularly  that  associated  with  mass  addition  normal  to  the  surface  of  injection. 
The  ultimate  objective  of  our  studies  is  to  demonstrate  that  the  long-time  (nonlinear) 
evolution  of  initially  small  acoustic  disturbances  is  intimately  coupled  to  the  specified 
boundary  forcing. 

Traditional  studies  of  acoustic  phenomena  in  models  of  solid  rocket  engines  (e.g., 
Culick  1990)  are  based  on  linear  and  weakly  nonlinear  stability  theory.  The  results 
of  this  work  have  been  shown  to  describe  many  of  the  events  observed  in  stable  and 
unstable  engines.  From  the  perspective  of  a  modeler,  the  formulation  decouples  the 
evolution  of  acoustic  amplitudes  from  the  driving  mechzmis.ns  associated  with  burning 
rate  transients.  It  is  also  notable  that  the  normal  mode  expansions  used  to  develop 
the  mathematical  model  are  valid  formally  only  for  a  totally  closed  chamber,  so  that 
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the  effect  of  the  nozzle  on  chamber  losses  is  not  properly  accounted  for. 

Our  work  emphcisizes  an  initial-boundary  value  approach  for  modeling.  An  ini¬ 
tially  steady,  low  Mach  number  flow  field  is  disturbed  by  specific  boundary  distur¬ 
bances.  Solutions  describe  the  totality  of  wave  systems  generated  by  the  disturbance 
and  their  evolution.  These  complete  solutions  are  in  contrast  to  those  based  on  the 
quasi-steady  cissumption.  In  particular,  we  can  describe  resonant  solutions  with  am¬ 
plitude  growth  that  are  not  part  of  the  solution  set  arising  from  the  quaisi-steady 
assumption.  It  is  notable  that  our  analyses  predict  the  appearance  of  relatively  com¬ 
plex  wave  systems  for  a  given,  more  elementary  boundary  disturbance. 

The  manuscripts  in  Appendices  A  and  B  describe  refraction  effects  on  the  acous¬ 
tic  time  scade  in  a  planar  channel.  On  this  time  scade  only  lineair  phenomena  are 
important  in  determining  the  wave  pattern  evolution,  a  conclusion  verified  by  the 
use  of  rational  approximation  methods  for  small  Mach  number  shear  flows.  The 
results  demonstrate  that  prescribed  axial  disturbances  are  altered  by  refraction  to 
create  higher  order  aociad,  oblique  and  even  purely  transverse  modes.  Hence  refraction 
provides  a  pathway  for  converting  energy  in  one  wave  mode  into  others,  depending 
on  geometrical  and  boundary  forcing  properties.  For  example,  an  aixial  boundary 
disturbance  in  a  semi-infinite  shear  flow  (traveling  waves)  can,  for  resonance  condi¬ 
tions,  produce  a  purely  transverse  wave  of  growing  amplitude.  This  type  of  result 
provides  a  source  of  understanding  about  the  appeau-ance  of  unexpected  wave  modes 
in  a  semi-confined  system.  Such  an  insight  could  not  be  obtained  from  more  tradi¬ 
tional  analyses  baised  on  the  quasi-steady  approximation  or  on  formulations  where 
only  one-dimensionaJ  (aocial)  processes  awe  assumed  to  exist. 

Work  in  progress  is  described  in  modest  detail  in  Chapters  2-4.  Acoustic  modes 
in  a  rectangle  with  one  open  end,  generated  by  sizable  wall  injection  transients,  are 
described  in  Chapter  2.  In  Chapter  3  we  consider  a  steady  injection- induced  low 
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Mach  number  flow  in  an  open-ended  tube  and  predict  acoustic  configurations  arising 
from  endwall  forcing  and  from  sidewall  injection  transients.  The  former  case  includes 
a  discussion  of  the  viscous  accommodation  layer  that  provides  a  smooth,  but  oscilla¬ 
tory  transition  between  the  acoustic  flow  and  the  no-slip  condition  at  the  wall.  An 
analog  problem  is  described  in  Chapter  4,  where  we  investigate  the  impact  of  the 
acoustic  properties  of  one  boundary  of  a  slab  when  the  other  boundary  exhibits  small 
amplitude  harmonic  motion.  The  eigenfunction  properties  are  shown  to  be  dependent 
on  the  explicit  conditions  chosen  and  that  is  shown  to  have  a  profound  effect  on  the 
evolution  of  the  entire  flow  in  the  semi-contained  gas. 

Engineering  design  of  reliable  solid  rocket  engines  can  benefit  from  an  awareness 
of  the  character  of  all  wave  systems  present  in  the  chamber.  The  primjtfily  analyticcd 
approach  presented  in  this  report  provides  a  sound  basis  for  prediction  of  wave  system 
behavior  on  the  acoustic  time  scale  of  the  engine  chaunber.  Additionally,  the  results 
provide  a  starting  point  for  the  consideration  of  longer  time  period  nonlinear  solutions 
that  describe  the  appearance  of  the  relatively  large  amplitude  response  of  low  Mach 
number  chamber  flows  to  smaller  prescribed  boundary  disturbances.  It  is  important 
to  recognize  that  typical  computational  solutions  to  model  rocket  engine  problems 
(e.g.,  Baum  1989)  do  not  provide  detailed  information  about  acoustic  wave  structure, 
but  rather,  give  only  point-wise  time  dependent  responses.  It  would  be  valuable  to 
develop  data  analysis  techniques  that  can  be  used  to  predict  acoustic  wave  structure 
from  nodal  point  data. 
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Chapter  2 


Acoustic  Wave  Modes  Generated 
by  Transient  Sidewall  Mass 
Addition 

2.1  Problem  Statement 

Sidewall  mass  awldition  is  used  to  mimic  the  gasification  of  a  burning  solid  propellant 
in  a  model  of  a  rocket  engine  chamber.  Velocity  injection  adong  the  sidewalls  of  an 
open  ended  rectangle  induces  a  low  Mach  number  internal  shear  flow.  A  positive 
transient  component  of  the  injected  mass,  superimposed  on  a  steauly  distribution  of 
the  same  magnitude,  is  the  source  of  initially  transverse  acoustic  wave  disturbances 
in  a  baisicaJly  steawly,  inviscid  rotationad  background  flow.  These  disturbances  are 
aissociated  with  an  aLxiad  pressure  gradient  which  acts  as  a  source  of  aociad  acoustic 
waves.  Reflection  processes  adone  in  the  semi-confined  open  ended  rectamgle  create 
distinct  aodad,  oblique  amd  transverse  wave  patterns.  Explicit  harmonic  mass  transient 
frequencies  aire  aissociated  with  rapid  amplification  of  specific  wave  mode  amplitudes, 
due  either  to  beat  phenomena  or  actuaJ  resonamce.  The  primary  objectives  of  this 
study  are  to  show  that  complex,  identifiable  wave  patterns  aue  created  by  relatively 
simple  boundary  disturbances,  and  that  boundary  driven  disturbances  can  be  the 
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source  of  laxge  acoustic  instabilities. 


2.2  Mathematical  Model 

Rational  perturbation  methods  are  used  to  derive  the  limiting  form  of  the  mathemat¬ 
ical  model  that  describes  acoustic  phenomena  in  a  low  Mach  number  {M  1),  large 
Reynolds  number  {Rt  »  1)  flow: 


p  =  pT, 

(2.1) 

pf¥  M  [(pu)j:  -f  {pv)y\  =  0, 

(2.2) 

p[ut  +  A/(uUr  -f  UUy)|  = 

(2.3) 

P  [ut  +  M{uvj;  -1-  uuy)]  = 

(2.4) 

pCy  [Tt  +  M(uTr  +  vTy)]  !=  -(7  -  l)Afp(u«  +  u»), 

(2.5) 

where  transport  effects  are  neglected  because  M/Rt  C  1.  The  axial  velocity  u  is 
nondimensionedized  with  respect  to  a  mean  induced  speed  =  v'^S  where  S  =  V  j  D' 
is  the  rectangle  aspect  ratio  and  the  characteristic  velocity  of  the  injected  fluid, 
is  used  to  nondimensionadize  the  transverse  velocity  component  v.  The  parameters 
M  amd  R^  are  defined  with  respect  to  U^.  The  axial  amd  transverse  dimensions  (z,  y) 
are  referred  to  (X',  W)  and  the  time  is  nondimensionalized  by  the  axiad  acoustic  time 
—  X'/Cfl  where  is  the  chaxaw:teristic  sound  speed. 

Mass  aiddition  on  the  sidewadls  is  described  by 

y  =  (0, 1):  V  =  (+, -)[uo(x)  +  K;(x,0]  (2.6) 

where  the  steady  and  tramsient  components  are  of  the  same  magnitude,  and  aire  both 
responsible  for  positive  mass  addition  only. 

The  induced  steady  internal  flow  arising  from  i;o(i)  adone  is  described  by  (p,  p,  T)  = 
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1  +  i\P{P,  RiO),  the  lowest  order  inviscid,  rotational  equation  system, 


UQx  -b  VOy  =  0, 

(2.7) 

1  „ 

UqUox  +  VoUoy  = 

- Pox, 

(2.8) 

7 

UoVqx  -f-  I'oVOy  = 

^  p 

- ^  0v> 

7 

(2.9) 

7  —  1 

eo  =  - — Po, 

II 

1 

O 

(2.10) 
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and  (2.6),  with  a  K,  =  0.  Eqs.  (2.10)  show  that  isentropic  conditions  prevail.  The 
wall  acts  as  a  source  of  vorticity,  so  that  the  no-slip  condition  uo(x,y  =  0, 1)  =  0  can 
be  satisfied.  Analytical  solutions  can  be  obtained  easily  for  ^  ^  1.  However,  the  flow 
and  pressure  fields  have  no  direct  effect  on  the  first  order  acoustic  waves,  ajid  are  not 
considered  further. 

2.3  Acoustic  Disturbances 

The  steady  flow  described  by  (2.7)-(2.10)  is  disturbed  by  turning  on  the  injection 
transient  ^^(x,^)  in  (2.6)  at  f  =  0.  The  mechanical  effect  on  the  boundary  creates 
a  velocity  field  (^,  V^)  and  thermodynamic  disturbances  defined  by  (p,  p,  T)  =  1-1- 
M{P,R,9).  One  may  reduce  the  equations  for  the  acoustic  fields  to  the  initial¬ 
boundary  vzdue  problem. 


Pou 

=  PQXX  +  S^PQyy, 

(2.11) 

t  = 

0;  Pq  =  Pot  —  0, 

(2.12) 

X  = 

0;  pQx  =  0,  X  =  1;  Po  =  0, 

(2.13) 

y  = 

0,1;  Poy  =  -j{-\-,-)Vy,t{x,t). 

(2.14) 

The  downstream  condition  in  (2.13)  is  chosen  for  simplicity.  A  more  physically  viable 
condition  might  be  Poi(l,  0  =  Q:^ot(l,  t)  where  a  may  be  thought  of  as  am  admittance 
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function.  The  complete  velocity  field  formula  given  by 


(u,u)  ~  (uo(x,y),uo(.r,y))  -  ^A)r,  j  df  (2.15) 

shows  how  the  wall  transient  enhances  the  steady  flow  field  described  by  (2.7)-(2.10). 
Additionally  9q  =  —  l)Ro  =  ((7  —  l)/'j)Po. 

Solutions  to  (2.11)-(2.14)  are  written  for  the  case  when 

V^{x,t)  =  ^-^[1  —  cos(u;t)]  cos  (2.16) 

7UJ  \  2  / 

which  describes  an  axially  distributed  non-negative  injection  velocity  with  a  harmonic 
time  variation. 


2.3.1  Nonamplified  Solutions 


For  u;  ^  n^o  s  [{^ ,  m  =  0, 1, 2,  •  •  •, 

Po  \ 

6^K  \2  J 

4wsin 


2  ^  cos(2mfl-y) 


[^(1)  —  u;2  ^mO 

u  cos  {2miry)  sin  (fimoO 


sin(u;t) 


(?)  - 


-“E 


UJ* 


mssl 


n„„(njL„-u>) 


(2.17) 


The  first  complete  term  within  the  parentheses  represents  the  quasi-steady  response 
of  the  gas  to  the  injection  transient.  The  second  term  can  be  written  cis  a  pair  of 
counter- propagating  purely  axiad  traveling  waves.  Finally,  e«M:h  of  the  components 
of  the  series  can  be  rewritten  as  a  pair  of  oblique  traveling  waves.  Clearly  one  will 
find  beats  when  u  —*  Q^no  ^d  rapid  amplitude  amplification  when  uj  =  fimOi  arising 
from  interaction  between  one  of  the  quasi-steady  modes  and  one  of  the  eigenvalue 
solutions.  These  possibilities  are  not  captured  in  a  traditionad  quasi-steady  analysis. 
This  elemetary  solution  displays  the  rich  complexity  of  wave  types  initiated  by  a 
relatively  simple  boundary  disturbance.  Even  more  complexity  will  appear  if  the 
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spatial  distribution  of  the  injection  transient  contains  more  axial  eigenfunctions  than 
used  in  (2.16). 

2.3.2  Amplified  Solutions 

W  hen  !jj  =  flmo,  m  =  0, 1,2,  •  ■  •,  (2.17)  fails  formally  and  one  must  seek  a  solution  in 
which  discrete  modes  are  amplified.  Fo*-  example,  when  ui  =  rr/2, 

j^  =  -|(cos(^()cos(^.)+...  (2.18) 

SO  that  the  amplitude  of  the  purely  axial  mode  (pairs  of  axial  traveling  waves)  grows 
linearly  with  time.  The  remaining  terms  in  the  solution,  represented  by  the  dots, 
describe  bounded  quasi-steady  and  oblique  modes. 

When  the  forcing  frequency  is  larger,  w  =  llmo,  nr  =  1,2,  •  •  •,  one  of  the  higher 
order  oblique  modes  is  amplified; 

=  -■^^cosiQrnot)  cos{2miry)  cos  y-ij  (2  19) 

where  the  dots  again  refer  to  bounded  modes.  Given  that  p  =  l-f-  MP,  the  amplified 
solutions  are  of  limited  V2due  on  the  time  scale  Mt  =  0(1),  which  happens  to  be  the 
characteristic  time  of  fluid  to  pass  through  the  semi-confined  container. 

2.3.3  Discussion  of  Results 

Table  2.1  (see  p.  12)  contains  results  for  a  system  where  =  L'/c^  =  lO'^s,  and 
6  =  5.  so  that  dimensional  frequencies  can  be  considered.  When  w'  =  125  Hz  the  re¬ 
sponse,  shown  in  Fig.  2.1  (see  pp.  13-16  for  figures  in  this  chapter)  for  Po/K  cos(7ri/2) 
IS  bounded  and  is  composed  almost  entirely  of  the  purely  axial  mode  and  the  y- 
independent  part  of  the  quasi-steady  solution  in  (2.17).  Beats  are  observed  in  Fig.  2.2 
when  a;'  s;  242  Hz,  the  response  £U'ising  primarily  from  the  same  modes  as  in  the  pre- 
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viously  discussed  case.  Linear  monotonic  amplitude  amplification  is  seen  in  Fig.  2.3 
for  the  purely  axial  mode  in  (2.18)  when  u;'  s;  250  Hz. 

The  bounded  result  in  Fig.  2.4  for  u'  sj  2500  Hz  at  y  =  0.1  is  primarily  due  to  the 
purely  a.xial  mode  in  (2.15)  but  includes  some  small  but  important  higher  frequency 
response  arising  from  y-dependent  parts  of  the  quasi-steady  part  of  the  solution. 

Quite  elongated  beats  are  seen  in  Fig.  2.5  for  u'  «  4992  Hz,  which  arise  from 
an  interaction  of  the  first  oblique  mode  with  the  second  mode  in  the  quasi-steady 
solution.  The  modulation  in  the  response  is  due  to  the  purely  axicd  wave. 

Finally,  when  u'  ^  5000  Hz  the  amplified  oblique  mode  in  (2.19)  is  found  in 
Fig.  2.6  at  y  =  0.1. 

The  results  demonstrate  that  relatively  simple  boundary  driven  disturbances  can 
produce  a  complex  acoustic  wave  response  in  a  semi-confined  geometry.  At  or  near 
the  resonant  driving  frequencies  explicit  modcd  amplification  occurs.  Even  in  the  case 
of  beats,  amplitudes  can  become  large  enough  during  the  growth  phase  to  render  the 
perturbation  approximation  invalid.  Clearly  one  must  resort  to  a  weakly  nonlinear 
theory  to  describe  acoustic  evolution  on  a  longer  time. 

Higher  order  linear  acoustic  theory  shows  that  in  addition  to  modal  amplification 
arising  from  boundary  driving,  other  sources  of  instability  include  quadratic  and 
cubic  acoustic  wave  interactions  and  refraction  of  the  basic  acoustics  by  the  shear 
flow.  The  latter  effect  creates  yet  more  wave  complexity  of  the  type  described  by 
Wang  cind  Kassoy  ( 1992a, b). 

The  results  described  here  aje  limited  by  the  simplified  boundary  condition  at  the 
open  end  of  the  geometry  in  (2.13)  and  by  the  lack  of  feedback  between  the  transien*^ 
wall  injection  and  the  pressure  field,  which  could  be  described  by  a  wall  admittance 
function  of  some  type.  The  pressure  node  boundary  condition  is  of  some  importance 
because  it  implies  that  no  work  is  done  at  the  exit  boundary.  As  a  result  the  acoustic 
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energy  evolution  in  the  cavity  depends  only  upon  the  work  done  on  the  injecaon 
surfaces. 

It  is  noted  that  the  acoustic  field  does  not  satisfy  the  no-slip  boundary  conditions 
on  the  sidewalls.  However  that  deficiency  can  be  easily  overcome  by  describing  a 
weakly  viscous  boundary  layer  adjacent  to  the  injection  surfaces.  This  layer,  resem¬ 
bling  an  injected  Stokes  boundary  layer,  is  significantly  thicker  than  a  traditional 
acoustic  boundary  layer  because  the  normal  velocity  arising  from  injection  is  "large" 
given  the  magnitude  of  U'^  (e.g.,  ^  Price  and  Flandro  (1991)  have  de¬ 

scribed  limited  aspects  of  this  boundary  layer  structure.  A  related  discussion  occurs 
in  Chapter  3  of  this  report. 

Each  of  the  limitations  noted  here  will  be  examined  in  future  endeavors.  In  partic¬ 
ular,  it  is  necessary  to  focuse  on  a  fully  interactive  formulation  of  the  model  problem 
which  includes  a  responsive  thermally  active  boundary  layer  (Kassoy  1979)  that  mim¬ 
ics  the  effects  of  the  combustion  of  solid  propellant. 


II 


Table  2.1;  .A.coustic  Response  Properties  for  Several  Forcing  Frequencies 


U! 

uj'  (Hz) 

Properties 

Primary  Response 

7r/4 

125 

stable 

axial  -r  y-independent  quasi-steady 
modes 

7r/2  -  0.05 

Auj'  ss  8  Hz 

beats 

(same) 

7r/2 

250 

axial 

amplification 

linear  growth 

2500 

stable 

axial  +  smaller  high  frequency 
modes 

flio  -  0.05 

Au)'  w  8  Hz 

beats 

first  oblique  mode  -b  quasi-steady 
mode  with  axi2d  wave  modulation 

5000 

oblique 

amplification 

linear  growth  of  first  oblique  mode 
-f-  eucial  wave  modulation 
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Fig.  2.6  The  pressure-time  response  for  w’»5000Hz. 
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Chapter  3 


Acoustic  Response  to  Boundary 
Disturbances  in  a  Self-Induced 
Low  Mach  Number  Shear  Flow  in 
a  Cylinder 

3.1  Problem  Statement 

Research  haa  been  conducted  to  study  acoustic  disturbances  in  a  low  Mach  number 
{M  C  1)  internad  flow  within  an  open  ended  cylindricaJ  chamber  of  length  L'  and 
diameter  D' .  The  mean  flow  is  generated  by  prescribed  side  wall  mass  addition  and 
the  disturbances  are  introduced  by  either  end  wall  transients  or  side  wall  injection 
transients.  The  objective  of  this  initiad-boundary  value  amalysis  is  to  find  the  com¬ 
plete  acoustic  response  arising  from  boundary  forcing,  including  all  possible  resonant 
modes. 


3.2  Mathematical  Model 

Perturbation  methods  are  used  to  derive  the  amalytic  form  of  the  mathematical  model. 
The  complete  set  of  fluid  mechanics  equations  is  nondimensionalized  with  respect 
to  the  characteristic  scales  of  the  variables  describing  the  flow  in  the  open  ended 
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cylindrical  chamber  with  low  Mach  number  (M  1)  and  large  Reynolds  number 
(Re  >  1). 


dp 

'idiprV,)  IdipVe) 

r  dr  ^  r  dB 

p 

'DVr 

V? 

1  =  f , 

r 

Dt 

r 

7  A/  dr 

'DVs  , 

v;v^si 

2  1  dP 

P 

.  Dt  ^ 

r 

^  -iMr  dB  ' 

dipK) 

dz 


=  0, 


_ _ 

^  dt  7A/  dz  ’ 

DT 

p_  =  -(,-l)A/P 
P  =  pT 


'lR(rK)  IdVg  d\\ 
r  dr  r  dO  dz 


(3.1) 

(3.2) 


(3.3) 


(3.4) 

(3.5) 

(3.6) 


where 


•dz)' 


The  transport  terms  are  neglected  since  (^)*^  <  1  for  most  of  the  situations.  The 
axial  velocity  component  is  nondimensionalized  with  respect  to  =  Vj,<5  where 
S  =  ^.  The  2izimutha]  and  reulial  velocities  are  nondimensionalized  with  respect  to 
which  is  the  characteristic  side  wall  injection  speed.  Time  t’  is  nondimensionalized 
with  respect  to  audal  acoustic  time  =  K  where  Cq  is  the  characteristic  sound  speed. 
Parameters  M  and  Re  are  defined  with  respect  to  The  mathematicad  model  can 
describe  a  combination  of  an  inviscid  rotational  steady  flow  with  a  similar  size  acoustic 
flow.  However,  the  description  of  viscous  phenomena  requires  more  general  equation. 


3.3  Steady  State  Flow 

.A  steady  state  mean  flow  is  generated  by  time  independent  mass  addition  at  the  side 
wall  that  mimics  propellant  burning.  In  the  parameter  range  of  (^)^Re  ^  1,  and 
\[  <C  1,  which  means  large  side  wall  injection,  all  the  variables  are  written  in  the 
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following  forms: 


P,  =  1  +  >PPo,  K  =  Ko,  I  =  1,2,3,  (3.7) 

with  P  representing  pressure,  density  and  temperature  respectively  and  K  represent¬ 
ing  radied,  azimuthaJ  and  axial  velocity  components  respectively.  The  leading  order 
equations  and  boundary  conditions  describe  an  inviscid  rotational  flow  that  satisfies 
the  no-slip  boundary  conditions  on  the  side  wall.  Analytical  solutions  for  the  radial 
and  axial  velocity,  as  well  as  the  pressure  distribution  for  ^  1  aire: 

Vor  =  sin(^r^),  Vq,  =  tt  /  K„(r)dr  cos(^r^),  (3.8) 

r  1  Jo  I 

Po  =  VmiT)dTVniz)dz,  (3.9) 

where  —Vn(z)  is  the  time-independent  side  wall  injection  distribution.  Related  solu¬ 
tions  can  be  found  in  Price  and  Fljmdro  (1991).  A  similar  iuaalytic  solution  can  be 
found  for  the  case  of  constant  side  wall  injection  when  ^  =  0(1).  The  flow  field  has 
to  be  obtained  numerically  for  any  other  form  of  the  side  wall  injection. 

3.4  Leading  Order  and  Higher  Order  Acoustics 

A  study  is  made  of  two  common  disturbance  sources  which  cause  an  unsteauiy  flow; 

1.  The  0{M)  end  wall  disturbamce  for  ^  »  1. 

2.  Transient  0{M)  side  wall  injection  for  ^  =  1. 

The  asymptotic  expansions  for  all  the  variables  ue; 

P  =  l-|-A/(Po  +  AfPi),  Vi  =  V,o  +  MVii,  t  =  1,2,3,  (3.10) 

where  P,  represent  thermodynamic  variables  and  K  represent  the  velocity  components 
in  r.d,z  direction  respectively. 
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3.4.1  End  Wall  Effects:  ^  1 

Leading  Order  Solution 


The  unsteady  part  of  the  leading  order  equations  can  be  combined  into  a  planar  wave 
equation  in  terms  of  the  axial  velocity  component: 


dt^  dz^ 


(3.11) 


The  corresponding  initial  and  boundary  conditions  are  given  by 


t  =  0, 

o 

II 

o 

II 

o 

(3.12) 

2  =  0, 

ViQ  =  sincjt; 

(3.13) 

2=1, 

dV^ 

(3.14) 

o 

II 

1- 

It  is  assumed  that  the  endwall  boundary  conditions  are  independent  of  r  and  9. 

-  t 

The  simplicity  of  the  equation  can  be  attributed  to  the  large  ^  ratio.  The  pre¬ 
scribed  velocity  at  2  =  0  hM  a  magnitude  similar  to  the  mean  ajcial  flow  velocity. 
The  boundary  condition  at  the  exit  is  derived  from  a  condition  of  zero  pressure  per¬ 
turbation  there. 

The  general  solution  for  Vto  is  : 


Ko(^,0  = 


sinu;t 


J...  {'I  [(A) 

-  {(^)  sia(A„.0  -f  <cos(6,..t)|sin(6„.r), 


(3.15) 


where  bn  =  {n  |)7r,  and  the  laat  term  represents  a  resonant  effect  present  only 
when  =  bn-.  A  careful  look  into  the  solution  provides  us  with  some  insight  into  the 
properties  of  the  acoustical  flow. 
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•  The  first  term  itself  and  the  second  part  of  the  nonresonant  Fourier  series  repre¬ 
sent  quasi-steady  motion  oscillating  at  the  driving  frequency.  The  other  Fourier 
series  terms  can  be  decomposed  into  two  counter-propagating  planar  travelling 
waves. 

•  If  the  driving  frequency  u  equals  none  of  the  natural  frequencies  the  solution 
is  bounded.  If  u)  is  very  close  to  one  of  those  natural  frequencies,  then  beat  will 
appear  due  to  the  interaction  between  the  quasi-steady  motion  and  one  pair  of 
travelling  waves. 

•  Resonance  occurs  if  u;  is  equal  to  one  of  the  and  the  amplitude  of  one  mode 
grows  linearly  with  time. 

Discussion  of  Results 

Table  3.1  contains  results  for  a  system  where  ^  =  10"^  s  so  that  dimensional 
frequencies  can  be  considered.  When  u>'  ss  159  Hz,  the  response  shown  in  Fig.  3.1 
(see  pp.  33-35  for  figures  in  this  chapter)  for  Ko  is  bounded,  and  the  contributions 
are  mainly  from  the  first  harmonic  term  and  the  first  few  Fourier  modes  in  (3.15). 
Beats  are  observed  in  Fig.  3.2  where  u'  as  238  Hz,  due  to  the  interaction  between  the 
quasi-steady  modes  and  the  first  axial  mode.  The  period  of  the  beat  is  approximately 
90  s.  The  linear  monotonic  amplitude  growth  seen  in  Fig  3.3  is  primairily  due  to  the 
resonant  axial  mode  in  (3.15)  when  u;’  =  250  Hz. 
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Table  3.1:  Acoustic  Response  Properties  for  Several  Driving  Frequencies 


u;'(Hz) 

Properties 

Primary  Response 

1 

159 

stable 

axial  -1-  quasi-steady  modes 

1.5 

238 

beats 

quasi-steady  modes 
with  axiad  wave  modulation 

7r/2 

250 

axial  amplification 

linear  growth 

Boundary  Layer  Correction 


The  leading  order  core  acoustic  solution  does  not  satisfy  the  no-slip  boundary  condi¬ 
tion.  A  viscous  layer  adjacent  to  the  walls  with  multiple  scale  structure  is  needed  to 
acconi  date  the  no-slip  boundary  condition  at  the  side  wall.  Experimental  data  (Price 
and  Flandro  1991)  implies  that  the  viscous  damping  thickness  is  Icirge  compared  to 
the  scale  of  spatial  oscillation  within  the  thin  layer.  A  multiple  scale  perturbation 
scheme  based  on  the  thickness  of  the  boundary  layer  and  the  characteristic  spatial 
length  of  the  oscillation  inside  the  boundary  layer  is  carried  out.  The  acoustic  solu¬ 
tion  in  the  core  provides  an  outer  boundary  condition  and  the  no-slip  condition  on 
the  side  wall  provides  an  inner  boundary  condition  . 

The  leading  order  acoustic  core  solution  shows  that  all  the  terms  can  be  classified 
into  the  following  two  forms:  V{z)e'^*  with  fl  =  u;  or  or  -t  cos(6„t)  sin(6n2)- 

The  viscous  layer  solution  for  the  outer  boundary  condition  =  V{z)e'^*  is  of 
the  following  form: 


r),  z,  t)  =  -V{z)  {  exp 


/ifl*  ifl 

.  ^  7T7~7Tv\  s 


(3.16) 


(V;.(^)) 

where  —  V^{z)  is  the  steady  side  wall  injection  velocity,  fi  is  the  dimensionless  viscosity. 
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The  variable  m  describing  spatial  oscillations  within  the  viscous  layer  and  n  describing 
viscous  layer  thickness  are  defined  as: 


where  a  =  M  and  3  = 
boundary  condition  V,q 


1  -  r 


a 


V  = 


3 


(3.17) 


Re{ 


The  solution  corresponding  to  resonant  outer 


=  —tccJhnt)sm{bnz)  is: 


=  -<sin(6n2)sin{fe„0 


+  {t  sin(^^)  sin{6„t)  +  bnt  cos(k()  cos(bnt) 


[  (  1  +  bn) 

V^z) 


2*^1 +6„) 


sin(A:^) 


sin(6„i) 


2{l+  bn) 
1 


^  cos{k^)  —  1  — 


^f>l 

1  bn 


siD(Ar^) 


cos(6„f)}  —  sin(6„z)e 


(3.18) 


where  k  =  ttVt. 

When  ^  =  r?  =  0,  the  solutions  satisfy  the  no-slip  boundaxy  condition  on  the  wall. 
On  the  other  hand,  when  ^  and  rj  — +  oo,  the  core  solution  is  recovered  in  an  oscillatory 
manner  since  the  amplitude  of  the  exponential  term  goes  to  zero  harmonically.  The 
effective  thickness  of  the  boundary  layer  depends  strongly  on  Q.  and  V^.  When  fJ  gets 
large,  this  thickness  will  become  small.  For  fl  =  u/,  if  u;  is  small,  the  viscous  effect  will 
penetrate  deep  into  the  core  region;  if  u;  is  large,  the  viscous  effect  will  be  confined 
in  a  very  thin  layer  near  the  side  wall.  The  same  effect  cam  be  expected  for  f2  = 
Thus  the  boundary  layer  thickness  is  smaller  for  higher  order  modes.  On  the  other 
hand  increasing  the  value  of  Vn,{z)  enhances  the  overall  boundary  layer  thickness. 

The  theory  is  valid  only  for  V'J/?  where  /?  <§;  1.  110  becomes  as  large  as 

of  0(1)  or  fl  is  small,  then  a  new  multiple  scale  perturbation  technique  is  needed  to 
find  solutions  where  weak  viscous  effects  fill  the  entire  chamber.  This  is  one  of  the 
future  research  topics. 
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Higher  Order  Acoustics 


The  Oi\P)  equations  can  be  rearranged  in  terms  of  Ki,Vni,Ki2  and  Pi,  where 
T’n  =  V;i{;,r.t),  V,i  =  ^11(2,  0  +  V;i2(2,  r,  t). 

5v;i2 


where 


dt 

IdirVn) 


=  9i{r,z,t), 

=  9'i{r,z,t), 


r  dr 

j  dt  ^  dz 

-tdz  dt  -^2(2, 
Pl  =  Px{z,t) 


gx{r,z,t)  =  K, 


dV, 


zQ 


dz 


2  dt 

f  f  \  ^  ^Po  \/  ^Po 


dz  dz 


/2(2,<)  =  - 


dt  ■  dz 

The  initial  and  boundary  conditions  are; 


dv^  -  a(K.,  +  Ko)  dv...  ,  „  dv,.. 

Po^r  +  *^*0 - - +  Vto— — t-  v; 


dr 


dz 


(3.19) 

(3.20) 

(3.21) 

(3.22) 

(3.23) 


(3.24) 


«  =  0, 

Uri  =  0;  Ki(z,r,t)  =  Pi(z,t)  =  0, 

(3.25) 

r  =  1, 

II 

0 

(3.26) 

2  =  0, 

II 

0 

(3.27) 

2=1, 

Pi=0. 

(3.28) 

In  addition,  K-iiVii  are  required  to  be  finite  at  r  =  0. 

Unlike  the  previous  lower  order  acoustic  solutions,  the  radial  and  axial  velocities 

for  this  order  depend  on  the  radial  variable  as  well  as  on  z  and  f. 

the  result  of 
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refraction  of  the  leading  order  acoustic  planar  waves  by  the  steady  state  mean  shear 
Howl  Wang  and  Kassoy  1992a).  The  driving  functions  ^i(r.z,t)  and  g2[r,z,t)  include 
the  products  of  mean  flow  terms  and  acoustic  terms,  and  quadratic  acoustic  wave 
interaction  terms.  The  latter  provide  a  source  for  acoustic  instabilities  in  addition 
to  boundary  driving  sources.  This  issue  will  be  addressed  and  discussed  in  the  later 
part  of  the  report. 

Comments 

This  study  of  the  acoustic  response  of  an  internal  flow  in  a  cylinder  to  0{M)  simple 
harmonic  velocity  oscillations  at  the  end  wall  reveals  that  significant  boundary  forcing, 
CIS  large  in  magnitude  as  the  basic  internal  flow  will  cause  numerous  travelling  wave 
modes  ais  well  as  a  quasi-steady  flow  motion.  Although  the  qucisi-steady  motion  is 
predominant  for  many  frequencies,  the  planar  travelling  waves  can  grow  significantly 
if  the  driving  frequency  is  close  enough  to  any  of  the  natural  frequencies.  Beats  can  be 
seen  under  those  conditions.  Resonance  will  occur  in  the  limiting  cases.  Higher  order 
acoustic  analyses  show  that  in  addition  to  the  boundary  forcing  sources,  interactions 
between  the  steady  mean  flow  and  leading  order  acoustic  waves,  and  interactions 
between  leading  order  acoustic  modes  provide  other  sources  for  acoustic  instability. 

When  resonance  occurs  or  beats  appear  for  some  caises,  the  linear  study  is  no 
longer  valid.  A  weakly  nonlinear  formulation  of  the  problem  is  needed  because  the 
perturbations  Me  not  small  for  sufficiently  long  time.  Detailed  discussion  is  deferred 
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to  later  sections. 


3.4.2  Transient  Side  Wall  Injection:  ^  =  0(1) 


Leading  Order  Solution 


In  this  case,  the  end  wall  velocity  is  zero  and  an  0{M)  transient  side  wall  injection  is 
superimposed  on  a  similar  magnitude  steady  state  side  wall  injection.  The  induced 
leading  order  acoustical  flow  field  in  the  core  is  described  by  the  following  multi¬ 
dimensional  wave  equation  and  initial  boundary  conditions  : 


d^Po  d^Po  1  dPo  1  d^Po  d^Po 

dt^  ~  drl  ^  ri  dri  r?  d9^  ^  dz^  ’ 


(3.29) 


i  =  0,  Pq  —  Pot  —  0, 

0  <  ^  <  27r,  Pq{6)  =  Po{9  +  2t) 

r  =  0,  jPol  <  oo, 

I  dPo  .f,  a 

«  =  0,  Poz  =  0, 

2=1,  Po  =  0. 


(3.30) 

(3.31) 

(3.32) 

(3.33) 

(3.34) 

(3.35) 


The  solution  for  a  general  positive  transient  side  wall  injection  distribution  hais  the 
following  form: 


Po{ri,9,z,t)  =  rx8f{z,9,t) 


OO  OO  OQ 


+  S  +  P,„„(t)sinm0]J„(u;„,„rx)cos(^^-^)7r2,  (3.36) 

/=1  msO  n=l  ^ 


where  ujmn  satisfy  =  0  and  rj  =  j.  The  time  dependent  coefficients  Aimn{t) 

and  Bimn{t)  are  determined  by  the  system: 


=  F,„„. 

t  =  0,  A,„.  = 


(3.37) 

(3.38) 
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"I"  ^^lmnBlmn{t)  ~  Gtmn, 


^  O'  Blmn  — *  ^S/mn(0)) 


dBlrnn 


5B(mn(0)t 


(3.40) 


where  —  (  2  A;mn>  Gtmni  5/t/mn(0),  5,4/mn(0)>  5fl/mn(0)  and  dgimni^) 

are  quantities  that  are  found  from  products  of  Fourier  series  expansions  of  each  de¬ 
pendent  variable. 

If  f{z,t)  =  cos(|2)  sin(u;t)  is  chosen  to  describe  an  axial  distributed  positive 
transient  component  of  the  side  wall  injection  which  is  superimposed  on  the  steady 
distribution  of  the  same  magnitude,  then  the  solutions  cire: 

•  Nonresonant  case;  u  ^  kn  for  all  positive  integer  n, 


Po  -  + 


E  •/o(wonri)  cos(|z)sin(u;0 


guait— steady  modes 
/lU)  ^  .  /’T  ,  /«■  , 

J— ^  sin(-0cos(-z) 


axtal  mode 


+  £  IT  ~  ITZ^l  sm(knt)cosi^z)Jo(u>onri).  (3.41) 


1=2 


oblique  modes 


The  first  term  represents  the  queisi-steady  response  of  the  gas  to  the  boundary 
forcing;  the  second  term  can  be  rewritten  as  a  pair  of  counter  propagating  purely 
axial  travelling  waves;  and  finally,  each  of  the  oblique  modes  can  be  rewritten 
cis  a  pair  of  oblique  travelling  waves. 


•  Resoneuat  case: 


=  f  (axial  resonant  modes), 


Po  =  -  —t  cos(^t)  cos{^z)  -I- 
W  Z  s 


(3.42) 


where  the  dots  represent  the  remaining  bounded  qucisi-steady  and  oblique 
modes. 

2.  ^  =  kri-,  kn  =  \/(f  +‘-0n 

Po  =  -{^tcos{oJt)Jo{!jJon‘ri)cos{^z)  +  . . . ,  (.3.43) 

2u;  I 

with  the  dots  representing  the  remaining  bounded  quasi-steady  axial  and 
oblique  modes. 

Discussion  of  Results 

Table  3.2  contains  results  for  a  system  where  =  ^  =  10~^  s  6  =  5  so  that 
dimensional  frequencies  can  be  considered.  When  J  ss  159  Hz,  the  response  shown 
in  Fig.  3.4  for  V^q  is  bounded  and  the  main  contributions  au'e  from  the  first  hcumonic 
oscillation  term,  the  first  few  ajcial  modes  and  quasi-steady  solution  in  (3.41).  Beats 
are  observed  in  Fig.  3.5  when  u)'  »  238  Hz.  The  period  of  the  beat  is  about  90  s 
due  to  the  interaction  between  the  quasi-steady  modes  cind  the  axial  mode.  Linear 
monotonic  amplitude  growth  is  seen  in  Fig.  3.6  for  the  first  axial  mode  in  (3.42)  when 
■jj'  =  250  Hz.  Higher  forcing  frequencies  will  activate  the  oblique  modes,  and  similar 
phenomena  as  just  discussed  can  be  observed  if  certain  range  of  forcing  frequencies 
are  picked. 

.\  more  detailed  discussion  of  related  solutions  for  acoustics  in  an  open  ended 
rectangle  can  be  found  in  Chapter  2. 
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Table  3.2:  Acoustic  Response  Properties  for  Several  Driving  Frequencies 


u; 

/'(Hz) 

Properties 

Primary  Response 

1 

159 

stable 

axial  +  qucisi-steady  modes 

1.5 

238 

beats 

quasi-steady  modes 
with  axial  wave  modulation 

7r/2 

250 

axial  amplification 

linear  growth 

Higher  Order  Solutions 

The  higher  order  analysis  is  done  to  discover  new  sources  of  2u:oustic  instability  above 


and  beyond  those  arising  from  the  boundary  driving  resonance  just  described.  The 

equation  system  for  the  higher  order  pressure  is: 

d^Pi 

dt^ 

,1  a/>,  1  ,  . 

dr\  ri  dri  rf  39^  dz^  ^ 

(3.44) 

f  =  0,  Pi  =  Pu  =  0, 

(3.45) 

0  <  ^  <  27r,  Pi{9)  =  Piie  -f  27r), 

(3.46) 

n  =  0,  |Pi|  <  oo, 

1  3P 

3  P 

z  =  0,  ~  = 

(3.47) 

(3.48) 

(3.49) 

z  =  1,  A  =  0. 

(3.50) 

There  are  two  nonzero  boundary  conditions  at  ri  =  ^  and  z  =  0  respectively. 
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The  solution  will  have  the  same  general  form  as  the  leading  order  acoustic  solution: 


Pi[rx,9,zJ)  =  X(ri,$,z,t) 

+  H  H  +  B/„„(<)sinm0]J„(u;,„„r)cos(^:4^)7r2  (3.51) 

^  =  l  m=:0  ri=l  “ 


where 


Ill  +  A'2  4I 

^  ‘/mn 

(0  = 

- 

*  /mn> 

t  =  0, 

- 

^Imn 

jAlmn 

(0), 

dt 

dB 

^  +  kUbU 

(0  = 

-  ^Ln. 

II 

o 

BLn  = 

SbIttiti 

iO), 

dBL^ 

dt 

II 

e 

.2/-1 

:  ( - IT 

^  2 

/  ‘  '^mn* 

^Almn(O); 


ffBlmn(0)i 


(3.52) 

(3.53) 

(3.54) 

(3.55) 

(3.56) 


^(mn.  ^i/mn(0).  ffBimr,(0)  and  gKimn(O)  ^Te  quantities  that  are  found 

from  products  of  Fourier  series  expansions  of  each  dependent  variable. 

The  forcing  functions  and  are  the  keys  to  the  solution  forms  of 
or  Resonant  growth  will  occur  if  or  Gi„„  contain  any  time  dependent 

harmonic  functions  which  have  the  same  frequencies  as  the  natural  frequencies.  Beats 
will  occur  if  forcing  frequencies  we  very  close  to  the  natural  frequencies. 


•  Nonreson2Lnt  leading  order  ui  ^  AOmn- 

The  forcing  functions  and  contain  (a)  numerous  terms  that  are  prod¬ 
ucts  of  steady  mean  flow  terms  and  the  leading  order  acoustic  terms.  The  latter 
consist  of  numerous  nonresonant  modes  having  the  same  frequencies  as  the  nat¬ 
ural  frequencies;  (b)  quadratic  products  of  the  leading  order  acoustic  terms  and 
terms  which  will  cause  0(t)  growth  for  1st  order  acoustic  modes  if  the  sum  or 
the  difference  of  the  two  leading  frequencies  or  the  boundajy  driving  frequency 
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and  the  leading  frequency  are  equal  to  one  of  the  natural  frequencies.  This  0[t) 
growth  is  not  from  boundary  forcing  directly. 

•  Resonant  leading  order  acoustics  =  Kimn- 

In  this  case,  the  leading  order  solution  contains  both  resonant  mode  and  non- 
resonaxit  modes.  Sources  for  the  time  dependent  growth  of  the  first  order  modes 
can  be  either  refraction  phenomena  or  quadratic  modal  interactions: 

1 .  The  highest  growth  rate  for  the  amplitude  of  the  first  order  acoustic  mode 
is  0{t^)  for  refraction  generated  disturbances. 

2.  For  quadratic  modes  the  contribution  to  the  largest  possible  growth  rate 
of  the  first  order  acoustics  can  be  divided  into  two  categories.  The  first  is 
from  the  product  of  one  nonresonant  mode  and  one  resonamt  mode  which 
leads  to  O(t^)  growth  if  the  sum  or  the  difference  of  the  resonant  frquency 
and  the  nonresonant  frequency  are  equal  to  any  of  the  natural  frequencies. 
The  second  is  from  the  product  of  two  resonant  acoustic  modes.  The  con¬ 
sequence  is  that  the  first  order  acoustic  modes  have  at  least  0{t^)  growth 
for  m2my  frequencies  or  even  O(t^)  growth  if  one  of  the  natural  frequencies 
is  twice  any  of  the  other  natural  frequencies  . 

The  conclusion  from  the  higher  order  analysis  provides  the  important  slow  time  scales 
of  different  conditions  so  that  weakly  nonlinear  analyses  can  be  carried  out  properly 
in  the  future  research. 

Comments 

The  analysis  of  the  acoustic  response  of  a  low  Mach  number,  injection  induced  flow 
ro  large  side  wall  mass  transients  has  been  done  for  the  inviscid  core  region  only.  .A 
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viscous  boundcury  layer  is  required  to  accommodate  the  no-slip  boundary  condition 
on  the  side  walls.  The  general  procedure  for  studying  the  viscous  layer  is  similar  to 
that  used  for  end  wall  case. 

.■\  weakly  nonlinear  study  of  the  acoustic  core  response  is  of  the  primary  interest 
because  the  results  will  predict  the  proper  long  time  behavior  of  the  disturbance 
evolution  and  these  nonlinear  core  results  will  provide  the  correct  outer  boundary 
conditions  needed  for  the  viscous  layer  analysis  on  the  longer  time  scale. 
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Chapter  4 

Effect  of  Boundary  Properties  on 
Resonant  Acoustic  Oscillations 

4.1  Motivation 

Traditional  studies  of  acoustic  instability  in  solid  rocket  engines  are  based  upon  the 
concept  of  normal  mode  expansion  (e.g.,  Culick  1990;  Price  and  Flandro  1991). 
Acoustic  wave  systems  are  represented  mathematically  in  terms  of  superposifons 
of  eigenfunctions  for  a  closed  chaunber  with  rigid  walls.  Such  a  representation  ig¬ 
nores  the  effect  of  finite  acoustic  admittance  of  the  burning  propell«uit  surface  and 
the  sonic  nozzle  entrance  on  the  acoustic  mode  shape.  The  applicability  of  normal 
mode  representation  to  rocket  combustion  chambers,  where  the  velocity  fluctuations 
perpendicular  to  the  bounding  surfaces  do  not  usually  vanish,  has  not  been  carefully 
examined. 

On  the  other  hand,  there  is  ^ul  extensive  literature  on  the  claissical  problem  of 
piston-driven  acoustic  oscillations  in  a  cylinder  (e.g.,  Ochmann  1981;  W^ulg  and  Kais- 
soy  1990,  Seymour  and  Mortell  1973a,b;  Jimenez  1973).  Results  show  that  at  res¬ 
onance.  the  nonline£u  limit  cycle  amplitude  and  wave  shape  are  strongly  dependent 
upon  the  imposed  boundairy  conditions  on  the  other  end  of  the  cylinder.  For  periodic 
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piston  motions  of  maximum  Mach  number  e  (e  <§;  1)  in  a  closed  cylinder,  shock  for¬ 
mation  occurs,  and  the  limit  cycle  amplitude  is  of  in  an  ideally  open  cylinder 

where  an  acoustic  pressure  node  is  present  at  the  end,  shock  waves  are  less  likely  to 
appear,  and  the  limit  cycle  amplitude  is  much  larger,  of  The  sharp  contrast 

suggests  that  boundary  conditions  have  a  profound  impact  on  the  behavior  of  an 
acoustic  system. 

It  is  the  purpose  of  this  study  to  explore  the  effects  of  different  boundary  conditions 
on  ga.sdynamic  wave  motions  in  a  confined  region,  and  to  develop  effective  mathe¬ 
matical  tools  capable  of  handling  this  type  of  problems.  A  simple  piston-cylinder 
configuration  is  used  as  a  paradigm  for  extracting  the  crucial  physical  phenomena  of 
interest.  Instead  of  the  traditional  quasi-steady  approach,  an  initiad-boundary  value 
approach  is  adopted  to  describe  the  time-evolution  of  the  amplifying  modes  to  limit 
cycles  under  resonant  boundary  conditions. 

4.2  Model  Problem  and  Formulation 

VVe  consider  the  one- dimensional  acoustic  wave  motion  in  a  cylinder  induced  by  an 
oscillating  piston  at  one  end.  The  condition  at  the  other  end  of  the  cylinder  is  specified 
by  an  acoustic  admittance  function,  or  response  function,  that  relates  the  disturbance 
acoustic  pressure  to  the  velocity.  The  p'ston  displacement  is  aissumed  small  relative 
to  the  cylinder  length,  and  the  characteristic  piston  Mach  number  c  «C  1. 

When  the  Reynolds  number  is  large,  the  flow  and  thermodynamic  properties  of 
the  cylinder  gas  are  governed  by  the  the  inviscid  Euler  equations  and  the  isentropic 
relations  (Wang  cind  Kassoy  1990): 


Pt  +  {pu)i  =  0. 

(4.1 

Ut  -1-  UUr  = 

(4.2 
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T  =  p^-\ 


(4.3) 


P  = 

The  thermodynamic  variables  p,  p,  and  T  in  (4.1)-(4.3)  are  nondimensionalized  with 
respect  to  their  equilibrium  values  po,  po  and  To;  the  velocity  u  is  nondimensionalized 
with  respect  to  the  equilibrium  sound  speed  Cq.  The  characteristic  length  and  time 
scales  used  in  the  above  equations  are  the  cylinder  length  U  and  the  acoustic  pessage 
time  ta  =  T'/coi  respectively. 

The  initial  and  boundary  conditions  are  given  by 

t=0,  p  =  p  =  T  =  l;  u  =  0,  (4.4) 

x  = - cos(u;t),  u  =  esin(u;t);  2:  =  1,  u  =  — (p— 1),  (4.5) 

ijj  'y 

where  the  dimensionless  piston  frequency  u  =  uj'U/c^.  In  the  second  boundary 
condition  the  acoustic  admittance  a  is  used  to  provide  a  genercd  relationship  between 
the  velocity  fluctuation  and  the  excess  (acoustic)  pressure  on  the  outflow  boundary. 
For  simplicity  we  assume  that  a  is  invariant  to  frequency  and  time.  One  notices  that 
(7  =  0  and  00  correspond  to  a  rigid  end  and  am  ideally  open  end,  respectively. 

4.3  Linear  Acoustic  Theory 

By  applying  the  perturbation  expamsions 

u  =  cui  +  c*U2  +  •  •  • ,  P  =  1  +  7(epi  +  H - ),  (4.6) 

to  (4.1)-(4.5),  the  following  describing  system  for  the  induced  leading  order  acoustic 
velocity  cam  be  derived: 

^itt  Uiix  ~  (4. 1 ) 

=  D)  =  Unit  =  0)  =  0, 

Ui(x  =  0)  =  sin(u;t);  uif(x  =  1)  =  -cruii(x  =  1). 
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(4.8) 

(4.9) 


This  system  is  most  easily  solved  by  using  the  Laplace  transformation  technique,  and 
the  solution  can  be  split  into  a  quasi-steady  part  and  a  transient  part; 


ui  —  -f-  UitTi  (■i-10) 

_  {sin  (w)  sin  [( 1  —  j)u;]  +  cos  (u)  cos  [(1  —  j)  cj]}  sin 
sin^  (w)  +  cos^  (a;) 

(Tsin(a>i)cos(u;t) 
sin^  (a;) -t- (T^  cos^  (u;)  ’ 

^ur  =  '3-^  — sin(Anx)e*^"L  (4.12) 

The  quasi-steady  solution  (4.11)  denotes  gas  oscillations  at  the  piston  driving  fre¬ 
quency,  and  is  the  long-time  solution  to  the  entire  system  provided  that  the  denom¬ 
inator  does  not  vanish.  The  transient  response  p^tr  given  by  (4.12)  is  a  summation 
over  all  the  eigenfunctions  multipled  by  time-dependent  coefficients.  The  eigenvalues 


An  are  defined  by  the  equation 


tan(A„)  =  icr. 


(4.13) 


They  are  complex  quantities  even  for  real  values  of  cr,  due  to  the  mixed  derivatives 
in  the  boundary  condition  (4.9).  By  writing  An  and  <7  as  An  =  An,  -t-  zA,;  cr  =  cr,  +  icr,, 
(4.12)  can  be  put  into  a  more  illuminating  form, 


(4.14) 


where 


An,  =  rZTT  -  0, 


A.  =  -iln  - - - 

2  1  +(7 


The  quantities  0n  and  0  are  both  functions  of  cr,  and  cr,,  whose  explicit  forms  are 
omitted  here  for  brevity.  Clearly,  (4.14)  represents  infinite  pairs  of  counterpropagating 
waves,  whose  amplitudes  attenuate  exponenti£dly  if  A,  is  positive. 

Several  observations  are  made  regarding  the  solutions  (4.10)-(4.15).  First,  for  real 
values  of  acoustic  admittance  cr,  the  imaginary  terms  of  the  solution  cancel  each  other 
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and  the  result  is  real.  Secondly,  special  result  for  hard  wall  boundary  (cr  =  0),  free 
boundary  (|cr|  — ►  oc)  and  infinitely  long  cylinder  (cr  =  1)  can  be  recovered;  in  the 
last  case  the  solution  reduces  to  a  simple  traveling  wave  of  the  form  sin[u>  (t  —  x)]. 
Thirdly,  it  can  be  shown  that  the  transient  wave  phenomena  decays  exponentially 
with  time  if  <7,  is  positive  and  less  than  infinity.  This  implies  that  the  acoustic  energy 
contained  in  the  initial  transients  are  rapidly  radiated  through  the  boundary,  and  only 
the  acoustic  mode  oscillating  at  the  piston  frequency  is  sustained  by  the  continuous 
mechanical  energy  imput  through  the  piston. 

Of  particular  interest  are  the  conditions  for  the  acoustic  system  to  become  res¬ 
onant.  This  happens  when  the  denominators  in  (4.11)  and  (4.12)  are  zero,  i.e., 
=  0  or  equivalently 

sin^(a;)  +  ^^cos^(w)  =  0.  (4.16) 

It  can  be  shown  that  (4.16)  has  solutions  only  if  (1)  <7  =  0  (rigid  wall);  (2)  |(7|  — >  cc 
(perfectly  open);  and  (3)  <7,  =  0  (<7  purely  imaginary).  When  cr  takes  other  values, 
there  exist  no  “right”  frequencies  to  make  the  system  resonamt. 

When  resonance  occurs,  the  linear  acoustic  theory  predicts  0(i)  growth  in  the 
leading  order  term  and  faster  growth  in  higher  order  terms.  The  solution  breaks 
down  after  some  time  elapse  and  nonlinear  theory  has  to  be  employed. 

4.4  Weakly  Nonlinear  Solutions 

.A.  systematic  perturbation  procedure  based  on  multiple-timescale  expansions  are  em¬ 
ployed  to  obtain  the  long-time  solutions  for  resonant  and  nonresonamt  acoustic  pro¬ 
cesses  in  the  cylinder.  This  work  is  currently  in  progress,  and  the  following  is  a  brief 
discussion  of  the  preliminary  results. 
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4.4.1  Resonance  in  a  Closed  Cylinder  (cr  =  0) 


The  linear  asymptotic  solution  at  the  resonant  frequency  u;  =  i/t  [u  =  integer)  reveals 
that 

u  =  ~et  [sin(u--j)  cos(u;t)]  +  e^t^f2{x,  t)  +  t)  +  ■  ■  ■ ,  (-I-IT) 


where  and  /a  are  bounded  functions.  The  asymptotic  expansion  breaks  down  when 
t  ~  u  ~ 

■A  two-timescale  expansion  based  on  the  fast  acoustic  time  t  and  the  slow  time 
r  =  are  employed  to  obtain  a  uniformly  valid  solution.  In  this  procedure,  the 
time  derivatives  in  the  original  equations  are  replaced  by  derivatives  with  respect  to 
both  t  and  r. 


X 


(4.18) 


and  u  and  p  are  expressed  asymptotically  as 


U  -f  eU'}  4- 


Solutions  for  the  leading  order  acoustic  quantities  ui  and  pi  consist  of  summations  of 
standing  wave  modes  with  slowly  varying  amplitudes  that  are  functions  of  r.  These 
amplitude  functions  have  to  be  determined  by  eliminating  resonant  terms  in  the  next 
order  solutions.  Under  null  initial  conditions  2ind  boundary  conditions  defined  by 
(4.5)  where  a  =  0,  we  obtain 


OO  30 

^  a„(r)cos(n7rt)sin(njrx),  p\  =  —  ^  an{r)sm{nirt)  cos{mrx}.  (4.20) 

n=l 

The  Fourier  coefficients  OnC''’)  are  calculated  from 


don 

dr 

don 

dr 


7  +  1 
16 


nx 


2  0!,fiQtn+m 


m=l 


E 


m=l 


^m^n— m 


-S 


nu^ 


(4.21) 


(4.22) 
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for  n  =  1,2,  The  quantity  is  the  Kronecker  delta  which  ensures  nontrivial 
solutions.  The  above  results  agree  with  the  solutions  obtained  by  Ochmann  (1985), 
who  used  the  method  of  averaging. 

The  infinite  mode-coupling  system  (4.21)  provides  a  mechanism  for  transfering 
acoustic  energy  from  the  resonant  mode  v  to  other  modes,  especially  the  high  fre¬ 
quency  modes.  As  a  result,  wavefronts  steepen  to  form  shock  waves  in  the  cylinder. 
The  internal  dissipation  of  the  shocks  cause  strong  damping  to  the  amplifying  wave 
field  and  eventually  a  limit  cycle  is  approached.  Numericcil  solutions  can  be  obtained 
by  truncating  the  infinitely  coupled  equations  (4.21)  to  a  desired  level  of  accuracy. 

4.4.2  Resonance  in  an  Open  Cylinder  (|(7|  — oo) 

The  eigenvalues  corresponding  to  an  open  cylinder  are  given  by  =  (n  —  1/2)t, 
n  =  1,2,  •  •  •.  The  i/th  mode  becomes  resonant  if  the  driving  frequency  w  =  A^.  In 
this  case,  based  on  the  behavior  of  the  linear  asymptotic  solutions  as  t  — ►  oo,  the 
correct  slow  timescale  needed  to  carried  out  multiple-scede  expansions  for  the  weakly 
nonlinear  wave  phenomena  is  r  =  and  the  limiting  amplitude  of  the  velocity 

and  acoustic  pressure  are  of  0(c‘^^).  The  larger  amplitude  arises  because  there  is  no 
shock  formation  and  therefore  less  energy  dissipation  prior  to  t  ~ 

The  proper  asymptotic  expansions  for  u  amd  p  on  the  time  scale  t  are  thus 

U  =  C3Ui  -f  CSU2  +  CU3 -I - ,  (4.23) 

p=  1 -|-7(c’Pi +c*P2  +  eP3  +  ---)-  (4-24) 

■Analogously  to  the  closed  cylinder  case,  the  leading  order  solutions  have  the  form 

OO 

“1  =  H  [a„(r)  co3(AnO  +  /3n(’')  sin(A„t)]  sin(A„i),  (4.25) 

'  n=l 

oo 

Pi  =  [-«r,(r)sin(A„f)  -f-/?n(r)cos(A„<)]cos(A„x),  (4.26) 

n=l 
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where  it  is  required  by  the  initial  conditions  that  a^(0)  =  J^(0)  =  0. 

The  slowly  varying  amplitude  functions  and  are  again  to  be  determined  from 
higher  order  considerations.  Unlike  the  usual  two-scale  expansion  processes,  where 
the  secular  equations  result  from  the  elimination  of  resonant  terms  in  the  next  order, 
solutions  to  are  found  to  be  bounded  and  thus  provide  no  clue  in  this  respect. 

Derivations  have  to  be  carried  out  to  0(e).  Eliminating  the  resonant  terms  at  this 
order  then  leads  to  the  requirement  that  a„  and  satisfy  an  infinitely  coupled  system 
of  first  order  differential  equations  like  (4.21).  However,  the  mode  coupling  terms  on 
the  right  hand  side  are  cubic  rather  than  quadratic,  because  quadratic  nonlinearity 
is  not  cumulative  when  a  pressure  node  is  present  at  i  =  1.  The  exact  forms  of  the 
amplitude  equations  are  currently  being  developed. 

4.4,3  Other  Considerations 

.•\s  discussed  in  Section  4.3,  another  type  of  acoustic  resonance  can  occur  when  a  is 
purely  imaginary  and  finite.  The  wave  amplification  and  nonlinearization  processes 
under  this  catagory  possess  new  features  that  are  also  worth  studying.  This  will  be 
undertaken  in  the  future. 

Also  of  interest  is  the  nonlinearization  of  nonresonant  acoustic  wave  systems. 
One  recognizes  that  the  three  types  of  resonamce- prone  boundary  conditions  men¬ 
tioned  earlier  aue  only  special  caises.  In  a  sense,  it  is  more  importamt  to  be  concerned 
with  neau-resonant  (beats)  amd  nonresonant  situations  as  well.  The  crucial  issue  is  to 
identify  the  conditions  under  which  shocks  form,  since  the  appearance  of  shock  waves 
introduces  am  additionad  damping  mechanism-internal  dissipation.  To  this  end  it  is 
imperative  to  use  the  correct  eigenfunctions,  given  the  boundary  acoustic  properties. 
Multi  pie- timescale  expansion  techniques  will  again  be  used  to  describe  the  wave  de¬ 
formation  process  due  to  the  nonlinear  accumulation,  which  occurs  on  a  longer  time 
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scale  in  comparison  with  the  resonant  cases. 


4.5  Concluding  Remarks 

An  initial-boundary  value  problem  is  formulated  to  study  the  linear  and  nonlinear 
evolution  of  acoustic  waves  in  a  cylinder  driven  by  a  small-amplitude  oscillatory  pis¬ 
ton.  The  study  serves  to  illustrate  the  drastic  influence  of  the  acoustic  properties  at 
the  cylinder  end  upon  the  time-evolution  and  limiting  amplitudes  of  the  confined  wave 
system.  The  mathematical  technique  developed,  based  on  the  multiple-timescale  per¬ 
turbation  expansion,  can  be  employed  in  the  future  to  investigate  similar  boundary 
effects  in  more  complicated  geometries  in  solid  rocket  engines. 

This  study  shows  that  the  acoustic  mode  shape  and  characteristics  are  very  sen¬ 
sitive  to  the  imposed  boundary  conditions.  Due  to  the  cummulative  nature  of  the 
quadratic  nonlinearity,  normal  mode  approximation  always  leads  to  the  formation  of 
shock  waves.  In  a  shock  free  system  this  generates  erroneous  results  and  overpre- 
dicts  damping  mechanisms.  More  exact  representations  of  the  system  eigenfunctions 
should  be  explored  in  the  future,  in  order  to  develop  reliable  prediction  capabilities 
for  rocket  engine  acoustic  stability. 
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Abstract 

A  systematic  perturbation  procedure,  based  on  small  mean  flow  Mach  num¬ 
ber  and  large  duct  Reynolds  number,  is  employed  to  formulate  and  solve  an 
initial- boundary  value  problem  for  acoustic  processes  in  a  shear  flow  contained 
within  a  rigid-walled  parallel  duct.  The  results  describe  the  general  transient 
evolution  of  acoustic  waves  driven  by  a  plane  source  located  at  a  given  duct 
cross-section.  Forced  bulk  oscillations  near  the  source  and  oblique  wave  gener¬ 
ation  are  shown  to  result  from  refraction  of  the  basic  planar  axial  disturbance 
by  the  shear  flow.  Refraction  also  causes  the  axial  waves  to  exhibit  higher  order 
amplitude  variations  in  the  transverse  direction.  As  the  source  frequency  ap¬ 
proaches  certain  critical  values,  specific  refraction  induced  oblique  waves  evolve 
into  amplifying  purely  transverse  waves.  .\s  a  result,  the  magnitude  of  the  re¬ 
fraction  effect  increases  with  time,  and  quasi-steady  solutions  do  not  exist.  The 
analysis  is  extended  to  the  thin  acoustic  boundary  layer  adjacent  to  the  solid 
walls  to  examine  the  shear  layer  structure  induced  by  the  variety  of  acoustic 
waves  in  the  core  flow.  Nonlinear  effects  and  acoustic  streaming  are  shown  to 
be  negligibly  small  on  a  scale  measured  by  a  few  axial  wavelengths. 


1  Introduction 


The  effect  of  shear  flow  on  acoustic  wave  propagation  was  first  studied  analytically 
by  Pridmore- Brown  (1958),  who  derived  the  following  linearized  wave  equation  for 
propagation  in  a  fully-developed  duct  flow; 

1  1  2.V/ 

~Pti  (1  A/  )prx  "b  Pyy  Pit  "b  2poCQMyVg,  (1) 

^^0  Co 

where  p  and  v  are  the  acoustic  pressure  and  normal  velocity,  respectively.  The  sound 
speed  of  the  mean  state  (po,  Po,  Tq)  is  Cq,  and  M  =  M{y)  is  the  shear  flow  Mach 
number.  Earlier  efforts  have  been  focused  on  seeking  quasi-steady  solutions  of  the 
type  p  =  Cross-stream  eigenfunctions  F  and  eigenvalues  k  are  obtained 

to  describe  the  shear  flow  distortion  of  specific  propagating  acoustic  wave  modes. 
Both  asymptotic  solutions  (Pridmore-Brown  1958)  and  numerical  solutions  (Mungur 
and  Gladwell  1969)  demonstrate  that  for  a  downstream  propagating  axi^d  wave  (the 
fundamented  mode),  the  acoustic  pressure  at  the  wall  is  significantly  larger  than  the 
value  at  the  centerline.  Calculations  for  upstream  propagation  (Hersh  and  Catton 
1971)  show  a  reversed  trend  of  acoustic  pressure  distribution. 

Quasi-steady  theory  is  useful  for  describing  only  limited  types  of  eu:oustic  phenom¬ 
ena  due  to  the  restrictive  nature  of  the  presumed  solution  form.  For  example,  one 
cannot  use  it  to  track  the  evolution  of  an  initial  disturbance  toweu^d  the  quasi-steady 
waveform,  if  it  exists.  Solutions  describing  temporal  amplitude  growth  (resonance) 
are  excluded  entirely.  Furthermore,  the  quasi-steady  solution  does  not  provide  the  ab¬ 
solute  magnitude  of  a  propagating  wave  and  its  relation  to  a  specific  acoustic  source, 
nor  does  it  include  new  waves  that  may  be  generated  by  refraction  of  the  given  wave. 
It  is  also  important  to  note  that  the  solution  as  well  u  the  formulation  exclude  the 
acoustic  boundary  layer  where  the  wave  motion  is  damped  by  viscous  effects  to  satisfy 
no-slip  conditions  on  the  duct  wall.  These  limitations  can  be  overcome  by  develop- 
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ing  an  initial-boundary  value  solution  for  acoustic  disturbances  in  a  shear  flow.  In 
addition  to  the  acoustic  analysis  one  must  consider  viscous  boundary  layer  effects 
adjacent  to  the  duct  walls. 

Previous  oscillatory  boundary  layer  analyses  are  mostly  for  incompressible  flows. 
For  example,  Stokes  (1851)  studied  the  long-time  quasi-steady  response  of  a  viscous 
fluid  to  boundary  oscillation;  Sexl  (1930)  and  Uchida  (1956)  investigated  laminar  pipe 
flow  due  to  oscillatory  pressure  gradient.  The  heat  transfer  process  in  the  pulsating 
pipe  flow  was  e.xamined  by  Romie  (1956).  These  studies  ali  demonstrate  the  velocity 
overshoot  at  the  edge  of  the  viscous  layer,  commonly  known  as  Richardson’s  annulax 
effect  (Richardson  and  Tyler  1929).  More  recently,  Barnett  (1970,  1981)  studied  the 
pulsating  pipe  flow  process  based  on  linearized  turbulent  Navier-Stokes  equations. 
Rott’s  (1980)  investigation  of  acoustic  oscillations  in  an  infinite  gas  region  parallel  to 
a  flat  plate  is  more  closely  related  to  the  present  study,  because  he  uses  a  low  Mach 
number  compressible  gas  model.  The  effect  of  mean  temperature  variation  along  the 
direction  of  oscillation  is  included,  but  no  mean  flow  is  allowed. 

In  contrast  to  the  traditional  quasi-steady  linear  approach,  Baum  and  Levine 
(1987)  developed  numerical  solutions  to  an  initial-boundary  value  problem  in  order 
to  describe  uni-directional  acoustic  propagation  in  an  axisymmetric  cylinder  with  a 
coexisting  mean  shear  flow.  The  code  is  based  on  Reynolds-averaged  Navier-Stokes 
equations  for  compressible  flow,  coupled  with  the  k-t  turbulence  model.  Acoustic  dis¬ 
turbances,  generated  by  a  disk-shaped  acoustic  source  of  spatially  uniform  strength, 
are  studied  over  a  few  acoustic  wavelengths.  In  this  short-time  calculation  one  cannot 
expect  to  find  the  quasi-steady  wave  structure  solution  used  to  solve  (1). 

The  present  study  is  inspired  by  the  limitations  of  the  classical  quasi-steady  solu¬ 
tions  mentioned  above,  amd  the  lack  of  long-time  results  in  Baum  and  Levine’s  work. 
The  physical  system  under  consideration  involves  a  horizontal  parallel  duct  contain- 
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ing  a  fully  developed  low  Mach  number  shear  flow  (cf.  Fig.  1),  a  two-dimensional 
duct  counterpart  to  the  cylinder  considered  by  Baum  and  Levine  (1987).  .An  initial¬ 
boundary  value  problem  is  formulated  for  an  acoustic  disturbance  propagating  into 
the  imposed  shear  flow.  The  disturbance  is  initiated  by  a  source  located  at  a  given 
duct  cross-section.  Such  an  approach  ensures  the  spontaneous  appearance  of  all  types 
of  acoustic  waves,  including  non-axial  waves,  arising  from  the  refraction  of  the  basic 
axial  wave,  and  provides  an  explicit  relationship  between  the  driving  acoustic  source 
and  the  evolving  wave  field. 

The  analysis  is  based  on  a  laminar  flow  model  for  a  viscous,  heat  conducting  fluid. 
By  using  a  rational  approximation  procedure,  in  Section  2,  transport  effects  are  shown 
in  a  formal  manner  to  be  limited  to  extremely  thin  acoustic  boundary  layers  adjacent 
to  the  duct  wall.  Perturbation  methods,  based  on  the  small  mean  flow  Mach  number 
parameter  .V/,  are  employed  to  find  solutions  for  both  the  transport- free  core  region 
in  Section  3  and  the  viscous  layer  in  Section  4.  The  solution  procedure  is  especially 
simplified  due  to  the  low  Mach  number  simplifications.  Finally,  in  Section  5,  the 
results  are  discussed  in  comparison  with  the  numerical  solutions  of  Baum  and  Levine 
(1987),  and  the  major  findings  of  the  present  work  are  summarized. 

The  results  demonstrate  that  refraction,  known  to  distort  the  pressure  distribution 
of  the  leading  order  axi^d  wave,  is  also  the  source  of  new  and  dispersive  acoustic 
transients.  When  nonresonant  conditions  prevail,  these  transients  evolve  into  oblique 
propagating  waves  and  a  forced  bulk  response  at  the  acoustic  source  frequency.  The 
former  correspond  to  the  selected  higher  modes  of  quasi-steady  propagation  in  duct 
acoustics,  while  the  bulk  response  is  composed  of  an  infinite  number  of  attenuated 
modes  that  decay  rapidly  away  from  the  plane  acoustic  source.  The  quasi-steady 
axial  wave  solution  agrees  with  those  from  the  classical  studies  (Pridmore- Brown 
1958;  Hersh  and  Catton  1971).  Resonance  occurs  in  one  of  the  propagated  modes  if 
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the  duct  width  is  some  integer  multiple  of  the  driving  acoustic  wavelength.  In  this 
case  purely  transverse  waves  with  growing  amplitude  are  found  to  exist  that  cause  a 
drastic  increaise  in  refraction  effects. 

In  the  viscous,  heat-conducting  acoustic  boundary  layers,  a  complex  response 
arises  from  the  variety  of  acoustic  waves  in  the  core.  The  transverse  velocity  in  the 
acoustic  boundary  layer  is  much  larger  and  more  complicated  than  that  predicted  by 
Rott  (1980),  because  it  must  match  with  the  core  solution  that  contains  refraction 
effects.  The  refraction  magnitude  and  the  acoustic  boundary  layer  thickness  obtained 
from  the  perturbation  solutions  are  found  to  be  comparable  with  those  of  Baum  and 
Levine  (1987). 

2  Mathematical  Formulation 


The  complete  dimensionless  equations  describing  the  compressible  fluid  motion  in  a 
planar  duct  shown  in  Fig.  1  can  be  written  in  the  form. 


p  =  pr,  (2) 

Pt  +  M  [(pu),  -f-  (pt;)v]  =  0,  (3) 

1  M  /  4  1 

P  +  iV/(uu*  -h  vUy)]  =  ~  \  3^  ^  3^ 

p  [u,  -f  A/(UU,  +  UUv)]  =  +  3“^ j  ’  (5) 


p[Ti  +  M{uTg  -h  v^y)]  —  — iV/(7  —  l)p(ui  -t-  Vy)  -H  'qp~^  ij^vy  "b  ^ 

(7-1)  [(«,  +  n'o.f  +  20'  (u'  +  -  jd'  (u,  +  »,)']  ,  (6) 

where  for  convenience  the  thermophysical  properties  are  assumed  constant.  The 
nondimensional  variables  are  defined  in  terms  of  dimensional  quantities  by 
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X  = 


f  ’ 


X  It' 

—  .  y  =  T- 

A‘ 


■) 


Quantities  pg,  Pg,  and  Tg  are  thermodynamic  properties  of  the  gas  at  mean  state,  d' 
is  the  width  of  the  duct  and  the  characteristic  velocity  of  the  mean  flow.  The 
characteristic  time  scale  is  defined  2ls  the  inverse  of  the  circular  frequency  of  the  axial 
acoustic  wave,  t'pi  =  1/uj',  so  that  the  wavelength  gives  the  proper  axial  length  scale 
=  d/ui'.  The  characteristic  transverse  velocity  v'^  =  Mui'd'.  The  dimensionless 
groups  in  (2)-(6)  defined  by 


.V/  =  ^,  R,  = 

d  y 


■'d' 


U  uJ'  ^ 

p.  =  -,  n  =  ^, 

a'  d 


(8) 


are  the  maximum  mean  flow  mach  number,  the  mean  flow  Reynolds  number,  the 
Prandtl  number,  and  the  normalized  axial  acoustic  frequency,  respectively.  One  no¬ 
tices  that  n  is  the  ratio  of  the  transverse  acoustic  time  in  the  duct  to  the  wave  period. 

In  the  present  study  solutions  to  the  above  described  system  are  sought  in  the 
limit  l/Pg  — *  0  and  A/  — ►  0,  where  it  is  aissumed  that  l/P*  A^-  Additionally,  one 
cissumes  that  the  Prandtl  numder  is  an  order  one  quantity  and  n  <  0(1). 

.As  in  previous  studies  (Pridmore- Brown  1958;  Mungur  and  Gladwell  1969;  Hersh 
and  Catton  1971),  the  basic  steady  flow  in  the  duct  is  assumed  to  be  fully-developed. 
It  is  driven  by  a  pressure  gradient  which  is  inversely  proportional  to  the  Reynolds 
number  of  the  flow.  One  can  easily  derive  from  (2)-(5)  that 


u  =  if{y),  u  =  0,  dpfdx  =  0{M^/ Re).  (9) 


Obviously,  on  the  length  scale  the  variation  of  p  is  negligibly  small. 

Since  the  Reynolds  number  is  very  large,  it  is  observed  from  (2)-(6)  that  except 
for  the  extremely  thin  acoustic  boundary  layers  adjacent  to  the  solid  surfaces,  which 
will  be  d’scussed  in  detail  in  Section  4,  the  wave  motion  in  the  core  region  is  basically 
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unaffected  by  transport  effects.  In  this  limiting  case,  the  state,  continuity  and  energy 
equations  (cf.  (2),  (3)  and  (6))  can  be  combined  to  give  the  familiar  results, 


These  isentropic  relations,  together  with  the  inviscid  version  of  (3)-(5),  suffice  to 
describe  the  acoustic  wave  motion  in  the  core. 

When  the  fully  developed  duct  flow  is  disturbed  by  an  0(c)  acoustic  velocity, 

u  =  U{y)->reu,  u  =  ct),  (11) 

it  can  be  shown  that  the  thermodynamic  corrections  are  always  O(Afc),  in  order  to 
balance  the  acoustic  components  in  the  governing  equations  (3)-(5).  Accordingly,  p. 


p  and  T  are  put  into  the  following  form: 

p=i  +  (a/c)7p,  p  =  i  +  (.v/c)p,  r  =  i  +  (A/c)r.  (12) 

The  continuity  equation  (3)  can  be  rewritten  in  terms  of  acoustic  variables  as 

^  -h  Us  +  Uy  +  MU{y)p^  +  [Me)  [(pu)r  +  (pv)y]  =0.  (13) 

Similarly,  the  i  and  y  momentum  equations  become 

ut  +  M  [U{y)u,  +  vU'iy)]  {Mt)iuur  +  vtiy)  =  -  ■  .  f .  (14) 

1  +  (Me)p 

vt  +  MUiy)v,  +  (A/c)(uus  +  vv,)  =  (15) 


In  (13)-(15)  terms  containing  U(y)  represent  the  shear  flow  interaction  with  the 
acoustic  field.  They  are  0{M)  quantities.  The  nonlinear  product  terms  are  of  0{eM). 
There  are  three  interesting  asymptotic  limits  that  can  be  applied  to  (13)-(15)  for  low 
Mach  number  shear  flows: 

I)  The  parameters  satisfy  the  inequality  e  .V/  1,  and  the  asymptotic  limit 

e  — ►  0,  ,Vf  fixed  is  used.  The  nonlinear  terms  can  be  ignored  since  their  magnitudes  are 
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small  relative  to  A/.  A  combination  of  (13),  (14)  and  (15)  generates  a  leading  order 
acoustic  equation  equivalent  to  (1).  This  demonstrates  that  the  previous  quasi-steady 
analyses  by  Pridmore- Brown  (1958),  Mungur  and  Gladwell  (1969),  and  Hersh  and 
Catton  (1971)  are  asymptotically  accurate  for  extremely  small  amplitude  acoustics 
relative  to  the  characteristic  shear  flow  Mach  number. 

2)  The  parameters  satisfy  the  inequality  0(e)  =  A/  ■C  1,  and  the  asymptotic  limit 
M  — »  0  is  used.  If  O(iVP)  terms  in  (13)-(15)  are  uniformly  ignored  relative  to  0{M) 
terms  in  the  aisymptotic  limit,  their  analogue  to  (1)  contains  no  0{M^)  term.  This 
shows  that  the  Pridmore- Brown  results  are  formally  valid  to  0{M)  where  M  <t.  1. 
.Numerical  studies  by  Baum  and  Levine  (1987)  are  primarily  concerned  with  this 
regime. 

3)  The  parameters  satisfy  the  inequality  M  e  <  0(1),  and  the  asymptotic 
limit  Sf  — *  0,  €  fixed  is  used.  An  examination  of  (13)-(15)  shows  that  the  nonlinear 
terms  and  the  refraction  producing  terms  containing  U(y)  are  of  the  same  magnitude, 
0(.V/),  in  the  limit.  As  a  result,  (1)  which  is  purely  linear  cannot  describe  acoustic 
phenomena  when  the  pressure  disturbance  in  (12)  is  0{M). 

The  following  analysis  is  focused  formally  on  Case  2  described  above,  although 
the  solutions  for  acoustic  variables  are  equally  v<ilid  for  Case  I  since  both  crises  yield 
identical  truncated  versions  of  (13)-(15)  to  0{M).  The  initial  state  of  the  disturbance 
quantities  is  described  by 


t  =  0,  u  =  V  =  p  =  0.  (16) 

.Acoustic  waves  are  excited  at  x  =  0  by  imposing  a  periodic  horizont2d  disturbance 
velocity 

X  =  0,  u  =  Asin(t),  (17) 


where  the  amplitude  .4  is  in  general  y-dependent,  although  the  simpler  case  of  con- 


stant  amplitude  is  emphasize  here.  The  normal  velocity  component  generated  in  the 
viscous  acoustic  boundary  layer  is  of  much  smaller  magnitude  relative  to  the  core 
flow  magnitude,  as  will  be  shown  later.  Hence,  so  far  as  the  acoustic  core  region  is 
concerned,  the  impermeable  condition 

y  =  0,l;  v  =  0.  (18) 

can  be  applied  directly  to  describe  the  acoustic  behavior  close  to  the  duct  walls. 


3  Core  Solution 

3.1  Acoustic  Transients  due  to  Axial  Wave  Refraction 

First,  the  acoustic  quantities  are  expanded  asymptotically  in  terms  of  M  in  the  fol¬ 
lowing  manner: 

-h  -f  =  (u,u,p,p,  r).  (19) 


In  order  to  account  for  the  small  changes  in  the  absolute  wave  propagation  speed  due 
to  the  0{M)  mean  flow,  it  is  necessary  to  use  a  strained  coordinate 


I  = 


1  -h  MUo  + 


(20) 


where 

Uo=rU{y)dy  (21) 

Jo 

is  the  bulk  shear  flow  velocity  area-averaged  across  the  duct.  One  finds  from  ( 13)-(  15) 
ordered  sets  of  equations, 


P\t  +  Uu  -h  V\y  =  0, 

(22) 

“u  +  P\i  =  0, 

(23) 

1 

Ui,-h  — piv  =0. 

(24) 
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and 


P2t  -T  U2z  -r  l'2i,  —  —i  {y)plr  l-Q<^lx- 


'  25 ) 


^2t  +  P2i  =  LoPli  -  i  iy)uii  -  i'Uyjvi, 
Q2‘ 


V2t  +  7riP2y  =  -f-’(y}Vli- 


Additionally,  the  isentropic  relations  (10)  imply  that 


26 1 


Pi  =  Pi. 


P2  —  P2- 


2S1 


Equations  (22)-(24)  are  combined  to  generate  the  linear,  homogeneous  wave  equa¬ 
tion  for  pi, 

Put  -  (piii  +  j^Pivv)  =  0-  (-9) 

The  leading  order  acoustic  equation  is  seen  to  be  unaffected  by  the  mean  shear  flow, 
except  for  the  bulk  convection  effect  incorporated  into  the  variable  x.  If  the  boundary 
velocity  oscillation  has  a  constant  amplitude.  .-4  =  1  in  (17),  the  solution  satisfying 
(16)-(18)  describes  a  wave  train  propagating  axially  into  a  quiescent  gas, 


Pi  =  =  sin(t  -  r).  n  =  0. 


(30) 


.Vote  that  ahead  of  the  wavefront,  when  i  >  t.  all  the  acoustic  quantities  are  zero. 

In  order  to  study  the  explicit  effect  of  the  shear  flow  velocity  gradient  on  the 
acoustic  field,  (25)-(27)  are  combined  to  give  the  second  order  analogue  of  (29), 

Pitt  —  (^Piit  +  ~  ~2f  oPir?  —  '-^'{y)p\it  +  ~^’^'iy)v\£.  (31) 

The  forcing  function,  representing  the  interactions  between  the  shear  flow  and  the 
leading  order  acoustics,  are  simplified  for  the  case  considered  here,  i.e..  when  the 
boundary  oscillation  is  y-independent.  Upon  inserting  (30)  into  (31),  the  latter  can 
be  rewritten  as 


Pitt  -  (  Pin  +  j^Plvv 


-2  \U{y)  -  f'oj  sin(f  -  x) 


I  TO  ; 
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for  X  <  t.  The  initial  and  boundary  conditions  necessary  to  solve  (32)  are  derived  by 
a  proper  combination  of  the  preceding  results  in  this  section  with  (16)-(18).  They 
are  given  by 


t  =  0,  P2  =  =  0,  (33) 

■r  =  Of  Pm  =  cos(f);  x -♦  oo,  pj  =  finite,  (34) 

y  =  0;  1.  pqj  =  0.  (35) 

The  hyperbolic  equation  system  {32)-(35)  can  be  solved  by  a  combination  of 
Laplace  transform  and  Fourier  series  method,  as  outlined  below.  If  Q  and  s  are  used 
to  denote  the  transformed  variables  of  pj  and  t,  respectively,  a  Laplace  transform  of 
(32)-(35)  with  respect  to  time  t  yields 

<3»  +  =  2  [C/(!()  -  t?„j  (36) 

x  =  0,  Qg  g  [L^(y)  -  U(i\  X oo,  Q  =  finite,  (37) 

y  =  0;l,  g,  =  0.  (38) 

The  homogeneous  boundary  condition  (38)  suggests  that  a  Fourier  series  solution  of 
the  form 

OO 

Q  =  ao(x,s)  +  ^  a„(i,3)cos(nTy)  (39) 

nsl 

is  obtainable.  Once  the  shear  flow  velocity  U{y)  is  Fourier  decomposed  into 

U{y)  =  UQ^-Y^UnCos{n-Ky),  (40) 

where  the  bulk  part  Uq  is  given  by  (21),  and 

f/n  =  2  /  f/(y)cos(nTy)dy,  (41) 

Jo 

it  becomes  clear  immediately  from  (36)-(38)  that  oq  =  0.  The  nth  Fourier  coefficient, 
for  n  =  1. 2,  •  •  is  governed  by 

-  (5^  +  <7^)an  =  (42) 

I  +  s' 
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where  the  parameter  is  defined  as 


HT 


The  solution  to  (42),  which  satisfies 


(43) 


where  Jo  is  the  zeroth  order  Bessel  function  of  the  first  kind. 

The  first  term  on  the  right  side  of  (46)  represents  a  quasi-steady,  axial  traveling 
wave  (fundamental  duct  mode),  with  y-dependent  amplitude  given  by  the  summa¬ 
tion.  It  is  an  0{M)  correction  to  the  leading  order  axial  wave  (cf.  (30))  and  thus 
describes  the  effect  of  shear  flow  induced  refraction  on  the  propagating  euciaJ  wave.  .As 
predicted  by  the  classical  theory,  the  axial  wave  acoustic  pressure  redistributes  itself 
nonuniformly  across  the  duct.  The  y-dependent  amplitude  function  grows  rapidly 
with  driving  acoustic  frequency  (cf.  (46)  and  (43))  and  varies  with  shear  flow  velocity 
profile.  It  is,  within  a  constant,  equivalent  to  the  result  of  Hersh  and  Catton’s  (1971) 
perturbation  study,  if  the  latter  is  rewritten  in  Fourier  decomposed  form.  However. 
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here  the  amplitude  is  completely  defined  because  the  solution  is  obtained  from  an 
initial-boundau'y  value  problem. 

A  second  fundamental  advantage  of  studying  linear  acoustic  refraction  phenomena 
in  terms  of  an  initial- boundary  value  problem  is  that  the  transient  (non  quasi-steady) 
evolution  of  the  acoustic  refraction  and  its  absolute  magnitude  can  be  found  explicitly, 
a  result  not  available  from  the  classical  studies  (Pridmore- Brown  1958;  Mungur  and 
Gladwell  1969;  Hersh  and  Catton  1971).  The  complete  second  term  in  (46)  represents 
acoustic  transients  initiated  by  the  passage  of  the  leading  order  axial  wave  through 
the  shear  flow,  including  dispersive  effects.  The  transients  evolve  into  oblique  waves 
(higher  propagated  modes)  and  forced  bulk  vibration  consisting  of  infinite  numbers 
of  attenuated  modes,  as  will  be  shown  soon. 

The  0(iV/)  axial  and  transverse  acoustic  velocities  can  be  obtained  by  integrating 
(26)  and  (27),  respectively.  They  are  listed  below: 


=  —  sin(f 


Uk  cosikiry)  +  ~  ^  j  ^  cos(nrj/), 

(47) 


(cos(l  -  i)  -  l| 


+  cos(C  -  { )  A  (?.\/{’  -  i‘)  sm(njrji).  (49) 


The  two  terms  in  both  uj  and  uj  expressions  again  represent  the  quasi-steady  axial 
wave  and  other  transient  phenomena,  corresponding  to  those  in  (46).  For  large  x. 
one  can  show  from  (48)  that  h„  C>(x"‘),  so  that  U2  remains  bounded  despite  the 
explicit  X  proportionality. 
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3.2  The  Evolution  to  Quasi-steady  Propagation 

Insightful  results  for  the  long-time  properties  of  pj  can  be  oucained  from  the  asymp¬ 
totic  properties  of  the  integral 


COs[t  -  0  d^. 


m 


VVhen  qn  ^  1.  (50)  converges  for  large  values  of  t  such  that 


lim  /n  =  < 

t—oo  I 


>  1 
qn  <  1 


(51) 


when  X  =  0(1)  (Gradshteyn  and  Ryzhik  1980).  It  follows  from  (46)  and  (51)  that 
the  nonresonant  long-time  solution  for  the  refractive  acoustic  pressure  can  be  written 

cLS 


»  2Uk 

P2  =  -sin(t  -  x)  — r-cos(iTy) 

fcTl  q  l 


-h 


f  a7  /2  s\n{t-^l-qlx) 

n=l  ''  Un  /  ^1  -  qi 


cos(mry) 


+  cos(0  On  — -l)— 7==-cos(nTy). 

n=.V+l  \^n  /  \Jq\-\ 


(5 


where  .V  is  defined  such  that  <  1  <  qs-\-\- 

The  second  full  term  in  (52)  contains  .V  Fourier  modes,  or  higher  propagated 
modes  in  classical  acoustics  terms.  Each  mode  can  be  rewritten  as  a  pair  of  oblique 
traveling  waves.  This  is  illustrated  by  rewriting  the  nth  mode,  denoted  by  as 


Pn  =  ^  [sin(t  -  Xni)  -I-  sin(t  -  X„2)]  , 


-  q 


(53) 


where 


'ni  =  y  1  -  X  -  HTy, 


-n2  =  )/ 1  -  9*  X  -h  riTry, 


(54) 
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represent  a  pair  of  oblique  paths  for  the  traveling  waves,  shown  in  Fig.  2.  The  wave 
reflects  repeatedly  from  both  duct  walls  as  it  travels  along.  Upon  each  reflection, 
it  switches  from  one  path  to  the  other.  If  the  transverse  coordinate  is  rescaled  by 
y  =  Qy,  so  that  both  V  and  j  are  normalized  using  the  same  characteristic  length 
scale  Xfl  =  (//>*;',  one  easily  finds  that  the  phase  speed  of  the  nth  mode  along  the 
paths  is  unity,  while  the  phase  speed  along  the  x-axis  varies  from  0  to  oo. 
depending  on  the  angle  of  incidence  =  sin“*  (^1  —  .  The  latter  is  identical  to 

the  angle  of  reflection  because  of  the  rigid-wall  assumption.  It  is  to  be  noted  that  when 
qn  is  close  to  1,  one  pair  of  large  amplitude  oblique  waves  become  nearly  transverse, 
so  that  a  form  of  wave  trapping  appears.  This  type  of  result,  to  be  discussed  in  full 
in  Section  3.3,  implies  that  a  resonance  occurs  when  qn  l~  and  amplitude  growth 
with  time  can  be  expected. 

The  oblique  traveling  waves  exist  only  when  qn  <  1.  In  dimensional  terms,  this 
implies  that  n  <  2d'/\'  (A'  is  the  wavelength)  must  be  satisfied  in  order  for  the  waves 
to  propagate.  The  number  of  non-axiai  traveling  wave  modes  is  thus  proportional 
to  the  duct  width  and  inversely  proportional  to  the  acoustic  wavelength.  This  is 
well  known  in  quiisi-steady  duct  acoustics  (see,  for  example,  .Morse  and  Ingard  1968). 
However,  the  present  transient  analysis  demonstrates  explicitly  that  refraction  of  a 
basically  aucial  wave  is  the  direct  source  of  the  oblique  propagating  waves.  These 
oblique  waves  will  also  interact  with  the  shear  flow  as  they  propagate  along,  to  gen¬ 
erate  more  complex  refraction  effects.  The  latter  are  not  included  in  the  solution 
because  they  are  0{M)  smaller. 

The  .ast  term  in  (52)  describes  a  bulk  response  of  the  gas,  driven  at  the  frequency 
of  the  acoustic  source.  These  so  called  attenuated  modes  {qn  >  1)  decay  exponentially 
along  the  x-axis,  and  thus  normally  affect  only  a  small  region  close  to  the  surface  of 
the  acoustic  source.  The  penetration  depth  is  proportional  to  (<7^  —  1)“U2  Given 
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(43),  one  finds  that  small  mode  number  n  and  high  frequency  Q  lead  to  deeper 
penetration.  Resonance  can  also  be  viewed  as  occurring  in  the  limit  (?„  — ♦  I"*",  so  that 
the  first  attenuated  mode  in  (52)  penetrates  asymptotically  far  into  the  field,  up  to 
the  bcisic  axial  wavefront.  In  practical  terms,  only  the  first  few  attenuated  modes 
from  the  second  infinite  summation  are  needed,  because  of  the  rapid  decay  of  the 
Fourier  coefficients  with  the  mode  number  n. 

The  transition  to  the  quasi-steady  solution  (52)  can  be  illustrated  by  numerically 
evaluating  the  second  order  acoustic  pressure  from  the  general  formula  (46).  Since 
the  Fourier  series  converges  fairly  rapidly,  only  the  first  20  terms  are  used  in  each 
summation.  comparison  of  the  results  with  those  from  summations  of  40  or  more 
terms  shows  agreement  to  within  three  decimal  places.  The  integral  /„  is  computed  by 
calling  the  QDAG  integration  subroutine  from  the  IMSL  software  library,  which  uses 
a  globally  adaptive  scheme  based  on  Gauss-Kronrod  rules.  Representative  examples 
of  results  for  various  acoustic  frequencies  and  different  types  of  duct  mean  flows  are 
discussed  below. 

Figs.  3  and  4  show  the  acoustic  refraction  effect  in  a  fully  developed  laminar  duct 
flow,  described  by  f/  =  4y(l  —  y).  In  Figs.  3a  and  4a,  the  time  variations  of  p^. 
evaluated  from  (46),  are  plotted  for  the  cases  of  fl  =  2  and  8  respectively,  on  a  duct 
cross-section  located  at  i  =  2t,  one  wavelength  downstream  from  the  plame  acoustic 
source.  The  solid  lines  represent  the  second  order  acoustic  pressures  at  the  wall, 
while  the  dashed  lines  denote  those  at  the  center-plane  of  the  duct.  For  comparison 
the  axial  wave  contribution  to  pj  (the  first  term  in  (46),  henceforth  denoted  as  p^j) 
corresponding  to  the  conditions  in  Figs.  3a  and  4a  are  depicted  in  Figs.  3b  and  4b. 
respectively.  The  mean  flow  Mach  r  umber  employed  in  the  calculations  is  0.1. 

Since  the  mean  flow  is  symmetric  with  respect  to  the  duct  center-plane,  the  Fourier 
coefficient  =  0  for  n  =  odd  (cf.  (41)).  Thus  the  first  oblique  wave  pair  for 
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p-2  corresponds  to  the  n  =  2  mode,  whose  cut-off  frequency  Q  =  (A'  =  d'  in 

dimensional  terms).  The  driving  acoustic  frequency  is  below  the  cut-off  frequency  for 
the  n  =  2  mode  in  Figs.  3,  so  that  only  the  propagating  axial  wave  exists  in  the  quasi¬ 
steady  state.  The  pj  curves  are  thus  similar  to  those  for  p^a-  The  small  discrepency, 
caused  by  the  acoustic  transients  induced  in  the  gas  medium  when  the  wavefront 
first  passes  the  given  location,  is  seen  to  die  out  gradually  as  the  solution  converges 
to  the  quasi-steady  solution  (52).  The  effect  of  attenuated  modes  is  negligibly  small 
at  X  =  2jr.  Calculations  conducted  at  other  frequencies  for  which  oblique  waves  are 
absent  show  that  the  transient  phenomenon  is  more  prominent  and  disappears  more 
slowly  for  lower  frequency  cases  than  for  higher  frequency  cases. 

The  situation  depicted  in  Figs.  4  is  quite  different.  Here  f)  =  8,  higher  than  the 
cut-off  frequency  for  the  n  =  2  mode.  The  second  order  acoustic  pressure  shown 
in  Fig.  4a  is  a  superposition  of  both  the  axial  wave  and  one  oblique  wave  pair,  in 
addition  to  the  small  transient  effect.  As  a  result,  it  is  dramatically  different,  both 
in  amplitude  and  phase,  from  the  pure  axial  wave  solution  oresented  in  Fig.  4b. 

Fig.  5  depicts  the  time  variation  of  pj  in  a  “turbulent”  mean  flow  field  emulated 
by  6'^  =  (1  -  |2y  -  The  other  plotting  conditions  are  identical  to  those  used 

in  producing  Fig.  3a,  the  laminar  flow  counterpart  of  Fig.  5.  The  refraction  induced 
acoustic  pressure  is  observed  to  be  much  smaller  in  amplitude  in  the  latter,  because 
the  mean  flow  represented  by  the  one-seventh  power  law  has  a  relatively  small  velocity 
gradient  in  most  parts  of  the  duct.  The  high  velocity  gradient  regions,  concentrated 
near  the  two  duct  walls,  are  too  narrow  to  promote  acoustic  refraction  on  a  global 
scale. 

The  bulk  convection  of  the  acoustic  wave  by  the  mean  flow  can  also  be  observed 
from  Figs.  3-5.  The  wavefront  emitted  from  the  plane  acoustic  source  arrives  at 
X  =  2ir  after  approximately  5.9  dimensionless  time  units,  which  is  less  than  2ir,  the 
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time  required  for  a  wave  to  travel  the  same  distance  in  a  static  medium. 

The  pressure  curves  in  Figs.  3  and  5  suggest  that  a  simple  relation  exists  for 
quasi-steady  acoustic  refraction  at  relatively  low  driving  frequency.  If  Q  <  2rr,  so 
that  non-axial  waves  are  absent  from  (52),  and  if  one  is  sufficiently  far  away  from 
the  plane  acoustic  source,  where  the  effect  of  the  attenuated  modes  is  negligible,  the 
entire  quasi-steady  solution  contains  axial  waves  only  (cf.  (19),  (30)  and  (52)).  In 
this  case  the  total  acoustic  pressure  normalized  by  its  value  at  the  duct  center-plane 
provides  a  quantitative  measure  of  the  global  refraction  effect: 

p,  ^  2i/k  fkir\ 

—  =  I  +  A/  —  -  cos(kTry)  ,  (55) 

Pic  fe=i  L  V  2  / 

where  the  subscripts  s  and  c  denote  quasi-steady  state  and  center-plane  respectively. 
Eqn.  (55)  is  in  fact  the  ratio  of  the  amplitude  functions  for  the  fundamental  mode, 
equivalent  to  those  found  in  earlier  studies  (Pridmore- Brown  1958;  Mungur  and  Glad- 
well  1969;  Hersh  and  Catton  1971)  in  the  low  frequency  range. 

In  Fig.  6  (55)  is  plotted  for  the  case  U  =  v  and  M  =  0.1,  for  three  types  of 
mean  flow  conditions:  U  =  1  —  12y  —  1|,  U  —  4y(l  —  y),  and  U  =  (1  —  |2y  —  1|)*^'. 
The  familiar  results  are  presented  in  decibels  to  follow  convention.  Obviously,  the 
acoustic  energy  of  the  downstream  propagating  axial  wave  train  is  channeled  towards 
the  walls.  The  linear  and  parabolic  mean  flows  cause  acoustic  refraction  effects  of 
similax  magnitudes,  while  the  shear  flow  represented  by  the  one  seventh  power  law 
generates  the  smallest  refraction  for  the  recison  explained  previously. 

It  must  be  emphasized  that  (55)  or  curves  like  those  in  tig.  6  are  accurate  rep¬ 
resentations  of  acoustic  refraction  phenomena  in  quasi-steady  state,  when  the  duct 
geometry  only  allows  purely  axial  wave  propagation  for  the  given  driving  acoustic 
frequency.  For  fl  >  27r,  like  the  case  depicted  in  Fig.  4,  refraction  of  the  planar  axial 
wave  also  generates  oblique  traveling  waves.  The  ratio  of  local  acoustic  pressure  to 
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that  at  the  center-plane  is  both  t  and  x  dependent,  and  is  not  a  useful  representation 
of  the  transverse  variation  of  the  refraction  effect. 

3.3  Resonant  Amplification  of  Refraction  Phenomena 

The  above  discussions  have  been  focused  on  nonresonant  situations.  The  refractive 
pressure  response  is  dramatically  different  if  =  i,  or  in  physical  terms,  d'  =  n’\'/2. 
where  n*  denotes  the  resonant  mode.  The  resonance  occurs  because  the  frequencies  of 
Jo(0  fhe  harmonic  function  in  (50)  are  nearly  identical  for  ^  »  1  when  q^.  =  [. 
.\s  a  result,  the  integral  becomes  unbounded  when  t  — ♦  oo,  and  no  quasi-steady 
solutions  exist.  When  the  time  is  large,  x  <  0(1),  one  can  show  that 

/„.(<>  l,x)  ~  cos  0  .  (56) 

In  particular,  at  x  =  0  (50)  can  be  evaluated  exactly  (Gradshteyn  and  Ryzhik  1980) 
to  give 

/n*(^0)=/  cos(t  ~  (5~) 

JQ 

which  agrees  with  (56)  if  the  aisymptotic  property  of  Jq  for  large  t  is  used.  By  using 
(56)  the  resonant  Fourier  mode  in  (46)  can  be  written  asymptotically  as 

Pn.{t  >  l,x)  ~  -bcos  -hfiy-  .  (58) 

It  shows  clearly  the  pair  of  purely  transverse  waves  trapped  in  the  duct  with  growing 
amplitudes. 

.Numerical  evaluations  of  (50)  for  qn  =  I  are  shown  in  Fig.  7,  where  the  horizontal 
coordinate  is  the  characteristic  coordinate  t  —  x  for  easier  comparison.  As  the  x  value 
is  increased,  similar  trends  of  growth  with  t  are  observed.  However,  the  solutions  at 
t  =  60  still  exhibit  strong  x-dependence  because  the  asymptotic  results  described  by 
(  56)  have  yet  to  be  reached. 


The  effect  of  resonance  on  acoustic  refraction  is  explicitly  demonstrated  in  Fig.  S. 
where  the  variations  of  pi  on  the  duct  wall  y  =  0.  at  three  separate  axial  locations  are 
plotted  against  time.  The  shear  flow  profile  is  parabolic  as  used  earlier.  The  driving 
frequency  Q  =  2t,  so  that  the  n  =  2  mode  in  the  second  term  of  (46)  is  resonant.  .All 
the  pressure  curves  show  growth  in  amplitude  after  sufficient  time  elapses,  when  the 
growing  resonant  mode  embedded  in  the  second  summation  dominates  the  system  in 
(46).  For  relatively  short  time,  however,  the  resonant  mode  amplitude  is  not  large 
compared  with  that  of  the  fundamental  mode  (the  entire  first  term  in  (46)).  The 
amplitude  of  pj  may  initially  decrease  with  time,  as  in  the  cases  for  f  =  2;r  and  3t, 
due  to  destructive  interference. 

Calculations  for  fl  values  slightly  below  and  above  the  resonant  frequency  2jr  show 
similar  trends  of  growth  for  p^.  However,  the  amplitudes  eventually  approach  large 
but  finite  limiting  values  predicted  by  (.52). 

The  resonant  result  illustrates  a  fascinating  mechanism  for  exciting  and  amplifying 
purely  transverse  waves  in  the  duct  through  axial  wave-shear  flow  interaction.  This 
shows  yet  another  distinct  advantage  of  the  present  analytical  based  initial-boundary 
value  study.  .Numerical  investigations  for  traveling  wave  refraction  (Baum  and  Levine 
1987)  are  limited  to  a  few  wave  cycles  only,  due  to  difficulties  associated  with  non- 
reflective  outflow  boundary  conditions.  The  results  in  Fig.  8  show  that  short-time 
solution  behavior  cannot  be  used  to  determine  if  resonant  amplifications  of  refractive 
pressure  are  occurring.  The  classical  qucisi-steady  solutions  (e.g.  Pridmore- Brown 
1958),  on  the  other  hand,  cannot  describe  the  resonant  mode  at  all. 

Due  to  the  growth  of  the  resonant  or  near-resonant  mode  .n  P2,  one  concludes 
that  for  such  systems  the  perturbation  expansion  (19)  breaks  down  as  t  ~  0(A/~^). 

new  derivation  will  be  needed  to  predict  its  long-time  behavior,  which  will  contain 
the  resonance-enhanced  refraction  effects  (the  transverse  waves)  in  the  leading  order. 
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It  should  also  be  noted  that  the  validity  of  the  above  developed  perturbation 
solutions  depends  on  the  value  of  0.  In  general,  once  the  dimensionless  frequency  fl 
becomes  as  large  as  the  magnitude  of  Mp^  becomes  comparable  with  that 

of  Pi  (cf.  (30),  (46)  and  (43)),  causing  the  breakdown  of  the  asymptotic  expansion 
(19).  The  actual  size  of  Mpj  nevertheless  also  depends  upon  the  magnitude  of  f/„. 
For  mean  duct  flows  symmetric  with  respect  to  the  center-plane  y  =  1/2,  the  largest 
terms  in  the  Fourier  summations  in  (46)  and  (52)  vanish  because  t/n  =  0  for  n  =  odd. 
As  shown  in  the  example  calculations,  pj  remains  0(1)  when  the  value  of  17  is  as  large 

3iS  8. 

In  the  high  frequency  limit  >  1,  when  acoustic  refraction  effects  are  no  longer 
small  correction  terms  of  0{M),  a  new  theory  needs  to  be  developed  which  includes 
refraction  in  the  leading  order  acoustics. 

3.4  Acoustic  Transients  due  to  y-dependent  Boundary  Dis¬ 
turbance 

In  this  section  non-axial  acoustic  transients  generated  by  source  oscillations  with 
y-dependent  amplitude  is  discussed. 

If  the  velocity  oscillation  at  the  source  (x  =  0)  is  given  by  u  =  i4(y)3in(t),  (29) 
must  be  solved  subject  to  the  following  conditions: 

<  =  0,  Pi  =  pi,  =  0,  (59) 

i  =  0,  Pit  =  -A{y)cos{t);  x  -►  oo,  pt  =  finite,  (60) 

y  =  0;  1,  Ply  =  0.  (61) 

A  solution  procedure  identical  to  that  used  to  solve  (32)-(35)  can  be  employed  to 
obtain 

JO 

Pi  =  Aoslait  -  i) +  AnIn{t,x)cos{niry),  x  <  t,  (62) 

nssl 
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where  is  the  same  integral  as  defined  in  (50).  The  Fourier  coefficients  .do  and  A., 
are  defined  by 

XI 

.4(y)  =  .do  +  .dnCOsfnTy).  (63) 

a=l 

Equation  (62)  consists  of  a  propagating  axial  wave  mode  of  constant  amplitude  .do. 
and  an  infinite  number  of  dispersive  higher  modes.  .\ote  that  unlike  the  0(.V/)  higher 
modes  due  to  refraction  considered  earlier,  here  the  higher  modes,  generated  directly 
by  the  acoustic  source,  are  of  the  same  order  of  magnitude  as  the  fundamental  mode. 
.After  sufficient  time  elapses,  in  the  absence  of  resonance,  the  solution  may  again  be 
expressed  into  the  qu«isi-steady  form, 


Pi 


=  ,dosin(<  -  x)  +  ^ 


.d 


•An 


sin  t 


+  cos{t)  Yl  - ^ 

n=.v+i  -  I 


qj^xj  cos(niry) 
e-\/?i-i^cos(njri/).  (64) 


.As  an  example,  the  boundary  disturbance  u(x  =  0)  =  jfsin(<)  in  a  laminar  mean 
flow  field  with  bulk  velocity  [/q  =  2/3,  is  considered.  This  bulk  velocity  corresponds 
to  U  =  Ay{l~y),  although  the  explicit  form  of  U{y)  is  not  needed  for  the  leading  order 
calculation.  The  maximum  mean  flow  Mach  number  is  assumed  to  be  0.1.  Figs.  9a-9d 
e.xemplify  the  characteristic  aw:oustic  pressure  signals  at  x  =  2t  (x'  =  A'),  under  four 
different  driving  acoustic  frequencies,  as  evaluated  numerictdiy  from  (62).  In  each 
figure,  the  time  variations  of  pi  on  both  duct  walls  {y  =  0, 1)  and  at  the  center-plane 
(y  =  1/2)  are  depicted.  In  Fig.  9a  fl  =  2,  lower  than  the  cut-off  frequency  for  the 
first  oblique  propagated  mode  in  the  duct.  The  resulting  wave  field  is  basically  an 
axial  one.  Thus  the  pressure  signals  at  three  different  y-locations  on  the  same  duct 
cross-section  are  almost  the  same,  whose  amplitudes  are  approximately  1/2  (.do  in 
(63)).  The  small  deviations  are  again  attributed  to  the  initial  transients  that  diminish 
as  time  progresses.  Fig.  9b  corresponds  to  fl  =  x.  the  resonant  frequency  for  the  first 
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acoustic  mode  (n  =  i),  which  causes  pi  to  grow  with  time.  The  wave  field  is  no 
longer  purely  axial  because  of  the  appearance  of  the  transverse  waves  associated  with 
the  resonant  mode.  The  results  in  Fig.  9c  for  Q  =  4  contain  both  the  fundamental 
mode  and  the  first  propagated  oblique  mode.  The  pressure  curves  vary  with  the 
y-coordinate  but  are  bounded.  Finally,  in  Fig.  9d  the  driving  acoustic  frequency  is 
high  enough  to  allow  the  third  propagated  oblique  mode  to  appear  in  the  duct,  in 
addition  to  the  fundamental  and  the  first  modes  (note  that  the  even  modes  do  not 
appear  because  A{y)  is  an  odd  function  of  y).  The  acoustic  pressure  curves  differ 
dramatically  from  those  in  the  previous  three  figures. 

Additional  numerical  evaluations  of  pi  show  that  when  the  acoustic  frequency  is 
varied  within  the  range  that  allows  a  fixed  number  of  propagated  modes,  the  character 
of  the  wave  field  remains  similar  though  the  results  vary  in  a  quantitative  sense.  How¬ 
ever,  whenever  the  cut-off  frequency  of  a  new  mode  is  crossed,  there  is  a  qualitative 
change  in  the  wave  phenomena. 

The  0(Af)  refraction  effect  of  this  more  general  acoustic  system  could  be  studied 
by  using  the  first  order  results  in  (31).  This  is  deferred  to  a  future  endeavor. 

4  Acoustic  Boundary  Layer  Solution 

4.1  Boundary  Layer  Formulation 

In  the  acoustic  boundary  layer  near  the  wall  at  y  =  0,  thermodynamic  perturbations 
must  be  of  the  same  order  of  magnitude  as  those  found  at  the  edge  of  the  core  flow. 
Therefore, 

p  =  1  +  Af^jp,  p  =  1  4-  A/^p,  r  =  1  +  Af^T,  (65) 
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where  the  tilded  variables  denote  acoustic  quantities  superimposed  on  the  mean  state, 
stretched  boundary  layer  coordinate  pointing  away  from  the  wall, 


n  = 


y 

6' 


(66) 


is  needed  in  order  to  describe  the  structure  of  the  extremely  slender  acoustic  boundary 
layer.  .A  balance  between  diffusion  and  other  important  physical  mechanisms  in  the 
general  governing  equations  (2)-(6)  can  be  obtained  if 


(67) 


which  provides  the  scale  of  the  acoustic  boundary  layer  thickness.  In  dimensional 
terms,  one  finds  the  well  known  result  S'  ~  indicating  explicitly  the  depen¬ 

dence  of  boundary  layer  thickness  upon  the  fluid  viscosity  and  the  frequency  of  the 
traveling  acoustic  waves.  The  proper  scaling  for  fluid  velocities  in  the  layer,  obtained 
by  examining  the  asymptotic  behavior  of  (11)  when  y  -♦  0  for  e  =  M,  are  given  by 


u  =  .V/u  +  (5r?6"(0)  +  •  •  • ,  u  =  (,V/(!»)v.  (68) 


The  0(8)  contribution  to  u  arises  from  the  Taylor  series  expansion  for  U{y)  in  the 
boundary  layer,  which  is  much  smaller  in  magnitude  than  the  horizontal  acoustic 
velocity  Mu.  given  the  assumption  made  in  Section  2  that  M  »  1/ Re-  Consequently, 
the  independent  variable  i,  rather  than  i,  is  an  appropriate  horizontal  coordinate. 

If  the  new  dependent  and  independent  variables  defined  in  (65)-(68)  are  used. 
(  2)-(6)  can  be  transformed  into 

p-t  =  M^pt,  (69) 

Pt  ^  Ur  +  f’n  =  [(Pfi)r  +  (pi>)nl  +  0{M6).  (  '0) 

U,  -  U,,,,  +  Pi  =  -M^ipUt  +  UUr  +  VU„)  +  0{M*)  +  0{M6).  M  n 

P.  =  0{6^). 


(73) 


+  (l'  -  '•)(Ur  T  I'r,)  =  {7(7  -  1)  -  p(!ir  +  7’t,) 

-{ptt  +  ar,  +  ot„)}  +  0(.v/^)  +  0{M6). 

The  above  equations  are  accurate  representations  of  (2)-(6)  in  the  acoustic  boundary 
layer  up  to  0(A/^).  The  orders  of  magnitude  of  terms  not  written  explicitly  are 
indicated  in  each  equation.  It  is  interesting  to  notice  that  nonlinear  product  terms, 
which  are  responsible  for  acoustic  streaming  phenomena  (Rott  1964),  again  turn  out 
to  be  of  O(M^)  relative  to  leading  order  acoustics.  Terms  of  0(\fS)  result  from 
the  residue  mean  flow  velocity  in  the  boundary  layer  (cf.  (68)).  These  are  higher- 
order  small  quantities,  of  comparable  size  with  0(M*)  terms  under  flow  conditions 
described  by,  for  example,  \f  =  0.1  and  Re  =  10®. 

In  the  following  solution  development  correction  terms  of  0(\P)  and  smaller 
in  (69)-(73)  are  ignored.  For  convenience  the  same  variable  names  (with  an  implicit 
zero  subscript)  will  be  used  to  describe  the  basic  acoustic  variations  in  the  asymptotic 
series.  .Notice  that  although  0{\f)  terms  do  not  appear  in  (69)-(73).  the  results  are 
valid  to  0(M),  and  are  matched  with  the  core  solutions  to  the  same  order.  The 
general  matching  conditions  are  expressed  mathematically  as 

-♦  00)  ~ 'Pcor<.(x,t,y  —  0),  'f  =  (u,  T).  (74) 

On  the  duct  wall  the  no-slip  condition  and  the  appropriate  thermal  conditions  must 
be  imposed. 

4.2  Transient  Solutions 

The  acoustic  pressure  is  seen  from  (72)  to  be  basically  uniform  across  the  boundary- 
layer.  equal  to  that  at  the  edge  of  the  layer.  Thus 


p  =  p(t.  x;  y  —  0). 


(7.5) 
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When  the  small  correction  terms  on  the  right  side  of  (71)  are  truncated,  the  resulting 
equation,  which  describes  the  transient  diffusion  process  of  the  horizontal  velocity 
perturbation  driven  by  the  acoustic  pressure,  is  seen  to  be  decoupled  from  the  others. 
It  must  be  solved  subject  to  the  no-slip  condition  at  r;  =  0.  At  <  =  0,  u  must 
vanish  because  there  is  no  acoustic  motion  in  the  core.  Following  a  standard  Laplace 
transform  procedure,  the  solution  is  obtained  in  integral  form; 

(76) 

Eqns.  (30),  (46)  and  (62)  can  be  used  in  (75)  and  (76)  to  evaluate  the  corresponding 
boundary  layer  velocity. 

It  is  also  of  interest  to  study  the  transverse  acoustic  velocity  and  other  thermody¬ 
namic  variables  in  the  acoustic  boundary  layer.  To  this  end  (69),  (70)  and  (73)  are 
combined  to  find 

ft  -  =  {j  -  l)pt,  (77) 

which  is  a  diffusion  equation  with  a  compressibility  forcing  function,  vaJid  to  0{M). 
An  appropriate  thermal  condition  needs  to  be  specified  at  the  duct  wall,  to  obtain 
solutions  for  T  as  well  as  p  and  v.  For  simplicity  two  idealized  types  of  thermal 
conditions,  i.e.,  the  adiabatic  condition  and  the  isothermad  condition,  are  considered, 
with  the  understanding  that  practical  situations  usually  lie  in-between. 

1)  Adiabatic  Wall  (T,(»/  =  0)  =  0).  A  formal  solution  to  (77)  subject  to  zero 
initial  condition  and  the  adiabatic  wall  condition  yields  the  isentropic  relations 

f  =  (7  -  l)p,  P  =  P  (78) 

if  (69)  is  used.  This  indicates  that  no  thermal  diffusion  exists,  and  all  the  thermody¬ 
namic  quantities  are  uniform  across  the  acoustic  boundary  layer.  Eqn.  (70)  can  be 
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integrated  to  0{M)  by  using  (76)  and  1 7S)  to  give 

0  =  -P,-,  +  I  <)«]  .  1791 

which  is  zero  on  the  wall  surface.  By  employing  integration  by  parts  for  the  error 
function,  the  above  result  is  rewritten  as 


=  -pt  + - ^ — =-  erf  i  — +  /  Pri  erf  ( -  tj 

+ - L_^  2 -  X  /  _  l')  +  ^  f  pTz\Jt-i  (e~^  - 

l  +  .V/^o  ^  V  /  ^  H 


where  the  factor  1/(1  +  MUq)  arises  from  dxidx.  Eqn.  (80)  is  arranged  according 
to  terms  that  grow  with  p  and  those  that  remain  0{l)  as  rj  oo.  The  former  are 
driven  by  the  refraction  induced  transverse  velocity  in  the  core,  and  can  be  shown 
to  match  with  the  asymptotic  behavior  of  the  corresponding  core  solution,  given  the 
coordinate  transformation  (66). 

2)  Isothermal  Wall  (T{r}  =  0)  =  Oj.  The  mathematical  system  for  T  is  analogeous 
to  that  for  u.  The  solution  can  immediately  be  written  down  as 

r  =  (7  -  1)^  P{(^,x)er/ (81) 

where 

f  (82) 

6t 

is  the  vertical  coordinate  for  the  temperature  boundary  layer,  of  thickness  character¬ 
ized  by  8t  =  [MH^PrRe))^-  Eqn.  (81)  can  be  transformed  into  a  more  meaningful 
form  by  integrating  the  right  side  by  parts,  and  defining 
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to  give 


^  =  (7-1) 


O  /"» 


p{t,x)  -  / 


(84) 


The  second  term  in  the  square  brackets  in  (84),  arising  from  conduction  effects, 
describe  the  deviation  in  temperature  from  the  isentropic  condition  represented  by 
p{t,i).  One  notices  that  at  i  =  f.  the  instantaneous  location  of  the  acoustic  wavefront 
in  the  core,  both  p  and  T  are  zero  to  match  those  quantities  in  the  undisturbed  flow 
field.  The  acoustic  density  and  transverse  velocity  for  the  present  case  can  be  derived 
from  (69)  and  (70),  respectively.  They  are  omitted  here  for  brevity. 


4.3  Quasi-steady  Solutions 


.As  shown  previously,  core  solutions  for  a  non-resonant  acoustic  system  as  t—x  becomes 
large  consist  of  quasi-steady  modes  only.  The  corresponding  boundary  layer  solutions 
are  more  easily  derived  by  using  the  complex  notation.  The  acoustic  pressure  is 
written  as  a  summation  of  complex  Fourier  series, 

p{t,x)  =  I'^One  (g5) 

n=0 

Given  (75),  a  comparison  of  (85)  with  the  core  pressure  expressions  shows  that 


no  =  l-.V/f 

*=i 

On  =  iW  - 


n  =  1,2,' 


for  the  axial  wave  refraction  case  (cf.  (30)  and  (52)),  and 


(86) 


ao  = 


On  = 


■Ao, 

•An 

V  I  ~9n 


n  =  1,2, 


(87) 


for  the  non-axial  wave  case  examined  in  Section  3.3  (cf.  (64)),  without  considering 
the  refraction  effect. 
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The  velocities  in  the  acoustic  boundary  layer  are  dependent  upon  the  boundary 
coordinate  r].  One  can  assume  that  the  quasi-steady  solution  for  each  horizontal 
velocity  mode  is  of  the  form  u„  =  6n(;?)c  ^  )  [f  this  expression  is  substitited 

into  (71)  with  (85),  and  the  conditions  Un(9  =  0)  =  0;  u„(r;  -+  oo)  =  0(1)  are  nvoked. 
then  the  desired  solution  is 


u  =  £  = 


n=0 


1  +  \fU, 


-  fl  (l  - 

0  n=0  ^ 


(88) 


where  On  is  given  by  (86)  or  (87),  depending  on  the  core  acoustic  solutions.  Eqn.  (88) 
can  be  rewritten  in  terms  of  its  real  part,  which  is,  for  the  refr<u:tion  case, 


u 


sin(t  —  x)  —  e 


=  ^_((,_,uf;H%) 
l+.Vt'oU  fl  / 

+  .V/j3^n(4“0  sin  (t-  yjl- 

nsl  \^n  /  L  ' 


-vf  f:  &,(|-.).-v^ 


n=.V+l 


sin(t)  —  e 


(89) 


As  the  value  of  i)  increases,  the  7-dependent  terms  in  (89)  diminish  exponentially, 
and  the  result  is  that  of  the  inviscid  core  evaluated  at  y  =  0. 

The  first  term  inside  the  curly  braces  in  (89)  is  associated  with  the  fundamental 
propagated  mode  in  the  core.  It  is  the  axial  traveling  wave  counterpart  of  the  classical 
Stokes  solution  (Stokes  1851).  Thus  the  characteristic  behavior  of  the  Stokes  solu¬ 
tion,  including  velocity  overshoot  (Richardson’s  annular  effect)  near  the  edge  of  the 
layer,  and  strong  viscous  damping  near  the  wall  will  be  observed  (cf.  Figs.  10).  The 
amplitude  of  the  fundamental  mode  deviates  by  0{M)  from  unity,  due  to  the  acoustic 
refraction  effect  (the  infinite  summation  in  the  parentheses)  and  the  bulk  convection 
(.V/f  0)  caused  by  the  mean  shear  flow  in  the  core.  In  the  case  of  downstream  propa¬ 
gation  the  former  effect  augments  the  axial  wave  amplitude,  while  the  latter  d<imps 
it.  For  upstream  wave  propagation  the  reversed  trend  is  obtained.  Eqn.  (89)  also 
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shows  that  the  .V  higher  propagated  modes  behave  in  the  same  way  as  the  funda¬ 
mental  mode,  except  for  decreased  amplitude  and  increased  phase  speed  as  the  mode 
number  n  increases.  Each  mode  exhibits  Rn.hardson's  annular  effect,  followed  by 
smooth  transition  to  no-slip  velocity  on  the  wall.  Like  the  core  solution,  the  effect  of 
attenuated  modes  are  limited  to  a  region  close  to  the  acoustic  source. 

The  result  given  by  (89)  is  illustrated  graphically  in  Figs.  10a  and  10b,  which 
depict  horizontal  velocity  profiles  across  the  acoustic  boundary  layer  when  Q  =  2 
and  7  respectively,  for  7r/4  intervals  over  one  acoustic  period,  at  a  location  one  and 
half  wavelengths  (x  =  3t)  downstream  the  plane  acoustic  source.  This  location 
is  sufficiently  far  from  the  acoustic  source  so  that  the  effects  of  stationary  modes 
(last  summation  term  in  (89))  are  virtually  non-existent.  The  solid  lines  denote  the 
horizontal  velocity  u  when  a  mean  flow  field,  described  by  f/  =  <.y(  l  —  y)  and  .V/  =  0.1, 
is  present  in  the  core  region.  The  velocity  distributions  for  the  case  of  no  mean  flow 
(.V/  =  0),  which  corresponds  to  the  Stokes  solution,  are  also  plotted  as  dashed  lines  for 
comparison.  It  is  observed  that  when  fl  =  2.  the  analogous  solid  and  dashed  curves 
differ  only  by  a  constant  multiplication  factor  (cf.  (89)),  because  the  only  propagated 
mode  in  the  duct  is  purely  axial.  The  velocity  amplitude  is  smaller  than  that  of 
the  Stokes  solution  due  to  the  bulk  convection,  whose  damping  effect  exceeds  the 
amplifying  effect  of  refraction  at  this  low  frequency.  It  should  be  pointed  out  that 
the  pairs  of  solid  and  dashed  curves  in  Figs.  10  are  plotted  at  the  same  rei  Hive  phase 
within  a  cycle  begining  at  their  respective  maximum  velocity.  The  two  solutions  are 
out  of  phase  in  the  absolute  sense  because  the  acoustic  wave  carried  by  the  mean  flow 
arrives  at  the  given  position  sooner  than  that  in  the  static  field. 

.As  the  driving  acoustic  frequency  becomes  higher,  the  increased  refraction  ampli¬ 
fies  the  velocity  oscillation.  When  fl  6.  the  refraction  effect  roughlv  balances  bulk 
convection  effect,  and  the  velocity  curves  are  found  to  coincide  with  those  derived 
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from  Stokes  solution.  The  acoustic  system  remains  dominated  by  the  single  funda¬ 
mental  mode  until  Q  >  27r,  wnen  the  second  mode  (n  =  2)  appears.  Fig.  10b  depicts 
such  a  case  where  H  =  7.  Here  the  amplitude  of  the  acoustic  velocity  is  larger  than 
that  of  the  Stokes  velocity.  The  velocity  profiles  can  no  longer  be  obtained  by  multi¬ 
plying  the  corresponding  Stokes  velocity  profiles  by  a  constant,  because  the  addition 
of  the  second  propagated  mode  alters  the  phase  of  the  velocity  at  each  time  instant. 
Additionally,  this  mode  makes  the  magnitudes  of  the  velocity  as  well  as  pressure  x- 
dependent  because  the  oblique  waves  strike  the  boundary  layer  nonuniformly  along 
its  course  of  propagation  (cf.  Fig.  2). 

On  the  acoustic  source  plaine,  x  =  0,  it  can  be  shown  that  (89)  takes  the  form 

u{x  =  0)  =  sin(f)  —  sin  - .  (90) 


\\1  the  0(Xf)  terms  disappear  because  acoustic  convection  and  refraction  take  place 
only  away  from  the  acoustic  source.  However,  (90)  does  not  satisfy  the  boundary 
condition  (17)  because  nonzero  v  is  allowed  at  x  =  0.  A  boundary-layer  type  of 
treatment  which  eliminates  the  slip  velocity  along  the  acoustic  source  plane  will  be 
necessary  in  order  for  the  extra  term  in  (90)  to  vanisn. 

The  acoustic  temperature  and  transverse  velocity  in  the  quasi-steady  state  de¬ 
pends  on  the  thermal  boundary  condition  along  the  duct  wall.  If  the  wall  is  adiabatic, 
the  thermodynamic  properties  of  the  boundary  layer  gas  again  obey  the  isentropic 
relation  (78),  and  the  transverse  acoustic  velocity  is  integrated  from  (70)  to  give 
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Figs.  11a  and  lib  display  the  transverse  velocity  profiles,  calculated  from  (91),  under 
conditions  identical  to  those  employed  in  Figs.  10a  and  10b  respectively.  The  trans- 
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verse  boundary  layer  motion  for  Q  =  7  is  quite  different  from  that  for  0  =  2.  owing 
to  the  existence  of  the  oblique  acoustic  waves  in  the  core  in  the  former  case.  .As  rj 
becom.es  large,  both  cases  exhibit  growth  in  velocity  amplitude  with  rj  to  match  with 
that  in  the  core,  in  contrast  to  the  transverse  velocity  in  a  Stokes  boundary  layer 
whose  amplitude  approaches  a  constant  as  the  edge  of  the  layer  is  approached. 

If  the  wall  is  kept  at  constant  temperature,  the  solution  for  t  and  v  can  be  derived 
from  (77)  and  (70)  in  the  same  manner.  They  are  expressed  below  in  the  complex 
form: 
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The  similar  forms  of  (92)  and  (88)  suggest  that  in  the  thermal  boundary  layer 
adjacent  to  an  isothermal  wall,  the  temperature  diffuses  in  the  same  way  as  the 
horizonal  velocity  does  in  the  viscous  boundary  layer.  The  first  two  terms  in  the 
curly  braces  of  (93)  are  identical  to  those  in  (91),  representing  the  effect  of  momentum 
diffusion  driven  by  the  acoustic  waves  in  the  core.  In  addition,  thermal  expansion  of 
the  gas  due  to  nonuniform  temperature  distribution  across  the  layer  also  contributes 
to  the  transverse  fluid  motion.  This  effect  is  represented  by  the  last  term  in  the  curly 
braces  of  (93). 

The  fundaunental  mode  in  (92)-(93)  can  be  shown  to  agree  with  the  results  ob¬ 
tained  by  Rott  (1980)  where  the  effect  of  the  mean  shear  flow  is  removed.  An  impor¬ 
tant  contribution  of  the  mean  sheaw  flow  is  to  generate  a  variety  of  acoustiv,  waves  that 
all  contribute  to  the  larger  transverse  velocity,  represented  in  the  acoustic  boundary 
layer  by  the  term  proportional  to  tj  in  (91)  and  (93).  In  the  limit  77  —♦  00,  the  mag- 
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nitude  of  this  term  exceeds  that  of  all  the  others,  thus  (91)  and  (93)  have  the  same 
asymptotic  behavior  which  can  be  shown  to  match  with  that  of  the  core  solution. 

5  Discussion  and  Conclusions 

In  this  study  a  systematic  analysis  has  been  developed  to  discover  the  effect  of  a 
low  Mach  number  shear  flow  on  acoustic  wave  propagation  in  a  planar  duct.  Two 
distinct  flow  regions  are  considered:  the  inviscid.  non-heat  conducting  core  region  and 
the  thin  acoustic  boundary  layer  near  the  wall  of  the  duct.  The  mathematical  analysis 
is  carried  out  in  the  limit  M  0  for  0(1/.V/),  when  the  axial  wavelength  is 

longer  than  or  comparable  in  the  order  of  magnitude  sense  with  the  duct  width. 

Solutions  for  the  acoustic  pressure  and  velocity  describe  both  short-time  acoustic 
transients  and  long-time  evolution  for  both  nonresonant  and  resonant  cases.  This 
study  bridges  the  recent  transient  numerical  study  of  Baum  and  Levine  (1987)  and 
earlier  qu£isi-steady  studies  (Pridmore- Brown  1958,  for  example)  and,  more  impor¬ 
tantly,  provides  new  results  not  available  in  those  investigations.  More  physical  in¬ 
sights  into  the  refraction  mechanism  are  obtained  by  demonstrating  explicitly  the 
interactions  between  the  mean  flow  and  the  various  types  of  acoustic  waves  repre¬ 
sented  by  Fourier  modes  whose  summation  describes  the  global  variations  in  acoustic 
quantities. 

It  is  of  interest  to  compare  the  present  perturbation  results  with  the  numerical 
solutions  of  Baum  and  Levine  (1987),  to  shed  light  on  a  number  of  issues  raised  in 
their  initial-boundary  V2jue  numerical  study  valid  over  a  few  acoustic  periods. 

1)  Acoustic  Refraction  Magnitude.  The  following  examples  are  used  to  demon¬ 
strate  that  the  present  linear  analysis  yields  refraction  magnitudes  comparable  with 
those  from  numerical  solutions  to  the  Navier-Stokes  equations  (Baum  and  Levine 
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1987).  In  a  duct  of  width  d'  =  0.1  m.  with  a  symmetric  mean  flow  described 
by  the  one-seventh  power  law,  a  center-plane  Mach  number  .\f  =  0.1  and  sound 
speed  d  =  340  m/s,  the  dimensionless  frequency  corresponding  to  /'  =  3000  Hz  is 
Q  =  2.772.  less  than  the  cut-off  frequency  for  the  first  oblique  wave.  The  quasi-steadv 
axial  wave  solution  is  thus  representative  of  the  acoustic  refraction  phenomena  since 
the  transient  effects  are  relatively  small.  The  pressure  amplitude  near  the  wall,  as  cal¬ 
culated  from  (55),  is  4.5%  larger  than  that  at  the  center-plane.  This  result  compares 
well  with  the  numerical  result  of  6.4%  by  Baum  and  Levine  (1987),  who  employed 
the  same  conditions  except  that  d'  is  the  diameter  of  an  infinite  cylinder.  When  the 
acoustic  frequency  is  changed  to  1000  Hz  (fl  =  0.924),  the  linear  asymptotic  solution 
and  nonlinear  numerical  solution  yield  near-wall  acoustic  pressure  increases  of  0.50% 
and  0.55%  respectively,  relative  to  the  centerline  pressures.  They  are  again  in  good 
agreement.  The  above  comparisons  should  of  course  be  interpreted  in  the  qualitative 
sense,  in  view  of  the  different  geometries  (parallel  duct  vs.  circular  cylinder)  and  flow 
models  (laminar  vs.  turbulent  k-t  model)  used  in  the  two  studies.  Nonetheless  they 
demonstrate  that  linear  studies  can  predict  refraction  effects  accurately,  and  that  the 
two  types  of  solutions  are  in  qualitative  agreement  if  comp£U’ison3  are  made  in  the 
same  piirameter  ramge.  In  the  light  of  these  conclusions,  it  is  likely  that  the  differences 
between  linear  and  nonlinear  results  noted  by  Baum  and  Levine  (1987)  result  from 
comparisons  in  inappropriate  parameter  regimes. 

2)  Acoustic  Boundary  Layer  Thickness  and  Structure.  The  boundary  layer  struc¬ 
ture  described  in  the  perturbation  solution,  including  the  Richairdson’s  annular  effect 
and  the  substantial  viscous  damping,  also  agrees  qualitatively  with  that  found  by 
Baum  and  Levine  (1987).  In  particular,  it  is  of  interest  to  compare  the  boundary 
layer  thickness  predicted  by  the  analysis  with  that  from  the  numerical  work. 

The  effective  thickness  of  the  velocity  boundary  layer,  as  defined  by  Lighthill 
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(1978),  is  given  by  56  —  5(.V//(n/?e))*^^.  According  to  Figs.  10,  this  corresponds  to 
the  distance  away  from  the  solid  wall  for  which  the  amplitude  of  u,  after  the  over¬ 
shoot,  approaches  its  asymptote  of  constant  value  to  within  approximately  2.7%.  For 
acoustic  waves  of  frequency  1000  Hz  traveling  in  air  contained  within  a  duct  of  width 
0.1  m,  under  standard  conditions  (one  atmospheric  pressure  and  room  temperature), 
the  calculated  boundary  layer  thickness  is  approximately  0.25%  the  width  of  the  duct. 
This  result  should  also  hold  for  wave  motion  in  a  circular  cylinder  of  diameter  equal 
to  the  duct  width,  because  the  curvature  effect  is  negligible  in  the  extremely  thin 
layer.  If  the  same  criterion  of  2.7%  deviation  is  applied  to  Fig.  14  of  Baum  and 
Levine  (1987),  one  finds  a  boundary  layer  thickness  of  approximately  0.35%  of  the 
diameter  of  the  cylinder.  This  result  is  actually  larger  than,  but  agrees  well  in  the 
order  of  magnitude  sense  with,  the  linear  perturbation  prediction.  Similar  agreement 
is  observed  in  terms  of  the  maximum  velocity  overshoots  and  the  locations  where 
they  occur.  Thus  we  do  not  agree  with  the  conclusion  of  Baum  and  Levine  (1987) 
that  linear  theory  significantly  overpredicts  the  boundary  layer  thickness. 

SJ  Nonlinear  Effect  and  Acoustic  Streaming.  In  the  present  work,  through  a  sys- 
temmatic  rational  approximation  and  perturbation  procedure,  it  has  been  demon¬ 
strated  in  (13)-(15)  that  the  convective  nonlinear  terms  are  0{M)  smaller  than  those 
responsible  for  acoustic  refraction  when  x  =  0(1).  Although  the  former  cam  have 
an  accumulative  effect  which  eventually  leads  to  waveform  deformation  and  weak 
shock  formation,  the  nonlinearization  process  becomes  prominent  only  adter  the  wave 
travels  a  distance  of  i  =  0(M~^)  (Kevorkian  and  Cole  1981;  Wang  and  Kassoy 
1990).  Nonlinearity  cannot  have  a  profound  influence  on  either  the  acoustic  or  the 
overall  flow  quantities,  on  the  0(1)  time  and  length  scales  considered  by  Baum  and 
Levine  (1987).  In  the  acoustic  boundary  layer  the  nonlinear  terms  are  again  shown 
to  be  0(.V/*),  relative  to  the  basic  veiriations  occurring  there  (cf.  (69)-(73)).  Acoustic 
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streaming  associated  with  the  nonlinear  convective  terms  is  thus  insignificant  relative 
to  the  amplitude  of  refraction  effects.  The  relatively  good  agreements  between  the 
present  fundamentally  linear  solutions  and  the  fully  nonlinear  numerical  soultions  in 
terms  of  refraction  size  and  acoustic  boundary  layer  structure,  discussed  above,  also 
attest  to  the  insignificance  of  the  nonlinear  phenomena. 

The  major  findings  of  the  present  study  can  be  summarized  £is  follows:  When  a 
plane  acoustic  source  of  uniform  strength  is  placed  across  a  duct  containing  an  undis- 
turbanced  shear  flow,  it  induces  leading  order  purely  axial  and  quasi-steady  acoustic 
waves  propagating  at  a  speed  which  is  modified  by  0{M)  due  to  bulk  convection. 
Second  order  acoustic  quantities,  including  j/-dependent  axial  and  oblique  propagat¬ 
ing  waves,  as  well  as  bulk  forced  oscillations  that  affect  only  a  narrow  region  near  the 
plane  acoustic  source,  are  generated  as  a  result  of  leading  order  axial  wave  refraction 
by  the  mean  flow  velocity  gradient.  The  propagated  and  attenuated  wave  modes  ex¬ 
hibit  transient  phenomena  initially,  and  evolve  gradually,  in  the  absence  of  resonance, 
into  their  respective  quasi-steady  state  long  after  the  passage  of  the  axial  wavefront. 

Resonance  occurs  when  the  duct  width  is  an  integer  multiple  of  the  the  driving 
acoustic  wavelength.  Then,  the  refraction  of  the  axial  wave  induces  an  amplifying 
purely  transverse  wave.  In  general,  the  refraction  effect  increases  with  the  driving 
frequency  as  well  as  the  mean  flow  Mach  number,  and  decreases  with  the  wave  ampli¬ 
tude.  The  refraction  induced  0{M)  wave  phenomena  become  increasingly  complex 
as  the  number  of  propagated  modes,  which  is  proportional  to  the  driving  frequency 
for  given  duct  geometry,  increases.  At  low  frequency,  when  the  driving  acoustic  wave¬ 
length  is  greater  than  the  duct  width,  the  only  propagating  waves  are  axial,  and  the 
net  effect  of  acoustic  refraction  is  to  distort  the  pressure  distribution  across  the  wave 
by  0[M).  The  quasi-steady  solutions  agree  with  the  classical  axial  wave  solutions. 

In  the  thin  acoustic  boundary  layer,  typically  with  thickness  of  less  than  1%  of 
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the  duct  width,  the  acoustic  pressure  is  basically  uniform  across  the  layer,  equal  to 
that  at  the  outer  edge  of  the  layer.  The  boundary  layer  responds  to  all  the  acoustic 
modes  existing  in  the  core  region,  generating  complex  velocity  and  temperacture  re¬ 
sponses.  Quasi-steady  solutions  again  exist  when  resonance  is  absent,  after  an  initial 
transient  period.  The  horizontal  velocity  component  for  each  acoustic  mode  exhibits 
Richardson’s  annular  effect,  followed  by  smooth  transition  to  no-slip  boundary  con¬ 
dition  on  the  wall.  The  total  horizontal  velocity  deviates  by  0{M)  from  the  Stokes 
solution,  because  of  the  acoustic  refraction  and  convection  effects  generated  in  the 
core.  The  transverse  velocity  grows  with  the  transverse  boundary  layer  coordinate, 
and  is  matched  by  the  core  solution  outside  the  layer. 
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Figure 

1.  The  physical  system:  source  generated  acoustic  oscillation  in  a  fully-developed 
duct  flow. 

1.  The  path  of  an  oblique  traveling  wave,  represented  by  (53)-(54). 

3.  Time  variations  of  (a)  the  second  order  acoustic  pressure  p^;  (b)  its  axial  wave 
contribution  pja  at  i  =  lir.  y  =  0  (duct-wall)  and  1/2  (center-plane),  as  calcu¬ 
lated  from  (46).  The  total  acoustic  pressure  p  =  Pi  -I-  .V/pj  for  the  wave  field 
generated  by  the  plane  source  disturbance  u(i  =  0)  =  sin(t)  in  the  laminar  duct 
flow  U  =  4y(l  -  y).  The  dimensionless  frequency  =  2,  and  the  maximum 
mean  flow  Mach  number  M  =  0.1. 

4.  Time  variations  of  (a)  the  second  order  acoustic  pressure  p?;  (b)  its  axial  wave 
contribution  p^a  it  x  =  2Tr,  y  =  0  and  1/2.  when  the  dimensionless  frequency 
n  =  8.  Other  conditions  are  identical  to  those  in  Figs.  3. 

0.  Time  variations  of  the  second  order  acoustic  pressure  p7  in  a  “turbulent”  mean 
flow  field  emulated  by  U  =  (1  -  i2y  -  The  locations  and  other  plotting 

conditions  are  identical  to  those  in  Fig.  3a. 

6.  .\xial  wave  quasi-steady  acoustic  pressure  profile  in  decibels  across  the  duct,  for 
three  types  of  shear  flows.  The  maximum  shear  flow  Mach  number  M  =  0.1. 
and  the  acoustic  frequency  fl  =  t. 

7.  Numerical  values  of  the  integral  (.50)  when  y,  =  1.  at  three  different  i-locations. 
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S.  Time  variations  of  the  second  order  acoustic  pressure  po  on  the  duct  wall  y  =  0. 
at  three  axial  locations,  at  resonant  frequency  H  =  2t.  Other  conditions  are 
identical  to  those  in  Figs.  3a. 

9.  Time  variations  of  the  first  order  acoustic  pressure  pi  at  three  j/-locations  on 
the  duct  cross-section  i  =  2x,  as  calculated  from  (62).  The  acoustic  waves 
are  generated  by  the  j/-dependent  velocity  perturbation  u(i  =  0)  =  ysin(t)  in 
a  mean  flow  field  described  hy  if  —  4j/(l  —  y),  \f  =  0.1.  The  dimensionless 
frequencies  are  (a)  fl  =  2;  (b)  fl  =  jt;  (c)  11  =  4;  and  (d)  11  =  10. 

10.  Profiles  of  quasi-steady  horizontal  velocity  across  the  acoustic  boundary  layer 
at  X  =  3jr,  for  ?r/4  intervals  over  one  acoustic  period,  when  the  mean  flow  in  the 
core  is  described  by  C/  =  4y(l  —  y)  and  A/  =  0.1,  amd  the  acoustic  disturbance  is 
u(x  =  0)  =  sin(<).  The  dashed  lines  correspond  to  the  Stokes  solution  (M  =  0). 
The  dimensionless  frequencies  are  (a)  fl  =  2;  and  (b)  11  =  7. 

11.  Profiles  of  quasi-steady  tramsverse  velocity  across  the  acoustic  boundary  layer 
adjacent  to  an  adiabatic  duct  wall,  at  x  =  3jr,  for  t/4  intervals  over  one  acous¬ 
tic  period,  under  conditions  identical  to  those  in  Figs.  10.  The  dimensionless 
frequencies  are  (a)  H  =  2;  and  (b)  11  =  7. 
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Abstract 

Acoustic-shear  flow  interactions  are  studied  in  a  rectangular  cavity  bounded 
by  impermeable  duct  walls  parallel  to  the  flow  direction  and  special  perpendic¬ 
ular  acoustic  reflectors  that  permit  the  passage  of  a  fully-developed,  low  Mach 
number  shear  flow.  Fourier  series-based  asymptotic  solutions  are  constructed 
to  provide  an  explicit  description  of  the  evolution  of  an  initially  imposed  axial 
standing  wave  disturbance.  The  bulk  convective  motion  of  the  shear  flow  is 
shown  to  be  responsible  for  periodic  axial  waveform  deformations.  .A.ddition- 
ally,  transverse  and  oblique  standing  acoustic  waves  as  well  as  single  frequency 
bulk  oscillations  arise  from  the  refraction  of  the  imposed  axial  acoustic  distur¬ 
bance  by  the  mean  flow  velocity  gradient.  Combinations  of  these  disturbances 
give  rise  to  irregular  acoustic  pressure  signals  in  the  duct.  Certain  refraction 
induced  transverse  and  oblique  acoustic  modes  are  amplified  under  resonant 
conditions.  It  is  shown  that,  in  the  parameter  ranges  of  solid  rocket  engines, 
the  refraction  and  bulk  convection  effects  are  small  in  general. 
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1  Introduction 


This  work  describes  the  spontaneous  evolution  of  a  spatially  distributed  small  initial 
disturbance  imposed  on  a  low  Mach  number  shear  flow  in  a  duct.  .An  initial-boundary 
value  formulation  is  used  to  predict  the  complete  spectrum  of  standing  wave  modes  (or 
equivalently,  combinations  of  traveling  waves)  generated  by  the  interaction  between 
the  initial  disturbance  and  the  shear  flow  in  a  finite  rectangul<ir  region.  The  primary 
objectives  of  this  study  are  to  demonstrate  that;  (1)  refraction  of  an  axial  wave  arising 
from  the  relaxation  of  the  initial  disturbance  is  the  source  of  a  myriad  of  higher  order 
(smaller)  transverse  and  oblique  acoustic  waves;  and  (2)  resonant  conditions  promote 
the  amplification  of  a  small  subset  of  refraction  induced  transverse  and  oblique  wave 
modes. 

Traditional  studies  of  acoustic-shear  flow  interaction  [1,  2,  3]  examine  the  quasi- 
steady  properties  of  a  sound  wave,  of  the  form  p  =  propagating  in  fully- 

developed  shear  flows  above  a  flat  surface  or  in  a  planar  duct.  Line£U'ization  techniques 
are  used  to  derive  the  fundamental  equation  governing  the  cross-stream  eigenfunction 
F  and  the  propagation  constant  k.  Asymptotic  solutions  [1,3]  and  numerical  solutions 
[2,  3]  predict  significant  distortion  of  specific  propagating  wave  modes  as  a  result  of 
acoustic  refraction.  When  the  wave  is  purely  axial  (the  fundamental  mode),  the 
refraction  effect  is  reflected  in  the  variation  of  F{y)  across  the  shear  flow  field. 

In  contrast  to  the  quasi-steady  linear  approach  used  in  Refs.  [1,  2,  3],  Baum  and 
Levine  [4]  used  numeric2d  methods  to  solve  an  initial- boundary  value  problem  based 
on  the  Reynolds-averaged  Navier-Stokes  equations  and  the  k-e  turbulence  model. 
Their  work  is  aimed  at  understanding  mechanisms  for  energy  exchange  between  the 
acoustic  and  mean  flow  fields  in  solid  propellant  rocket  engines.  In  the  model  problem 
discussed  in  Ref.  [4],  an  initially  steady  shew  flow  field  in  a  rigid  walled  axisymmetric 
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cylinder  is  disturbed  continuously  at  an  inlet  or  outlet  plane  to  generate  unidirectional 
traveling  acoustic  waves.  The  waves  move  only  a  few  acoustic  wavelengths  before 
computations  are  terminated,  so  that  the  qutisi-steady  solution  in  Refs.  [1.  2,  3]  is 
not  attained  on  the  short  time  scale  involved. 

More  recently,  Wang  and  Kassoy  [5]  used  an  initial-boundary  value  approach  to 
consider  the  two-dimensional  planar  duct  counterpart  of  the  refraction  problem  stud¬ 
ied  in  Ref.  [4].  A  systematic  perturbation  procedure  based  on  the  small  mean  flow 
Mach  number  \f  is  used  to  demonstrate  that  nonlinear  effects  are  0{M)  smaller 
than  the  refractive  contribution  to  the  total  acoustic  pressure.  The  Fourier  series- 
based  solution  in  Ref.  [5]  describes  short-time  acoustic  tramsients  arising  from  axial 
wave  interaction  with  the  shear  flow,  as  well  as  their  evolution  into  long-time  quasi¬ 
steady  forms.  In  the  nonresonant  case  the  solution  includes  not  only  the  usual  axial 
propagating  mode,  but  also  a  finite  number  of  oblique  propagating  modes,  and  an  in¬ 
finite  number  of  non-propagating  bulk  modes  that  decay  rapidly  away  from  the  plane 
acoustic  source.  When  resonamce  is  present,  a  pair  of  amplified,  trapped  transverse 
waves  also  appear.  These  distinct  modes  have  not  been  extracted  from  the  numerical 
computations  in  Ref.  [4]. 

It  is  of  interest  to  extend  the  aforementioned  traveling  wave  models  of  acoustic- 
shear  flow  interaction  to  axially  confined  geometries  in  which  multiple  wave  reflections 
occur.  This  is  done  in  the  present  work  by  considering  standing  waves  trapped  be¬ 
tween  two  wave  reflectors  separated  by  a  finite  axial  distance  in  the  duct  examined  in 
Ref.  [5].  The  idealized  model  is  used  as  a  paradigm  to  demonstrate  the  surprisingly 
complex  response  of  an  initially  steady  shear  flow  to  an  imposed  axiad  velocity  dis¬ 
turbance  in  the  confined  region,  and  to  infer  the  acoustic  convection  and  refraction 
magnitudes  in  solid  rocket  engine  chambers.  The  inclusion  of  aocial  wave  reflection 
effects  in  this  study  provides  improved  relevance  to  rocket  engine  acoustics,  in  com- 
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parison  with  the  previous  traveling  wave  studies. 

Fourier  series-based  analytical  solutions  give  explicit  modal  response  that  may 
not  be  ecisily  extracted  from  numerical  data.  The  initial  disturbance  generates  a 
leading  order  acoustic  field  which  is  purely  axial.  It  interacts  with  the  shear  flow  to 
generate  0(A/)  acoustic  disturbances  consisting  of  purely  axial  standing  waves,  tran- 
verse  standing  waves,  oblique  standing  waves,  and  bulk  fluid  oscillations  (Helmholtz 
modes).  The  linear  combination  of  these  disturbances  gives  surprisingly  irregular 
0{\I)  pressure  signals  at  a  given  location  that  could  be  mistakenly  identified  as  a 
chaotic  or  turbulent  response.  The  ratio  of  initial  disturbamce  wavelength  to  duct 
width  and  the  relationship  between  the  disturbance  and  resonant  frequencies  are 
found  to  be  important  factors  affecting  the  refraction  magnitude,  which  is  in  general 
small  except  under  resonant  conditions.  When  resonance  is  present,  amplifying  modes 
exist,  which  can  be  represented  as  pairs  of  transverse  and/or  oblique  traveling  waves 
with  growing  amplitudes.  The  results  of  the  present  analysis  are  employed  to  show 
that,  in  the  parameter  ranges  of  typical  solid  rocket  engine  gasdynamics,  refraction 
effects  can  produce  only  0(10%)  variations  in  acoustic  pressure. 

2  Mathematical  Formulation 

The  conceptual  model  of  the  present  acoustic-shear  flow  interaction  study  is  schemat¬ 
ically  illustrated  in  Fig.  1.  A  fully-developed  shear  flow  sweeps  through  a  p«irallel  duct 
region  of  length  L',  confined  by  rigid,  impermeable  side  walls  at  y'  =  ±d’.  The  ther¬ 
modynamic  state  of  the  steady  flow  field  is  defined  by  (p(„  pj,,  Jq),  and  the  equilibrium 
speed  of  sound  c(,  =  ^'tRTo,  where  R  is  the  gas  constant  and  7  the  ratio  of  the  specific 
heats. 

.\n  axially  distributed,  y'-independent  initial  disturbance  to  the  steady  flow  field 
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evolves  into  transient  acoustic  oscillations  on  the  acoustic  timescale  =  L' /c'^.  which 
is  much  shorter  than  the  mean  flow  passage  time  because  the  maximum  mean  flow 
Mach  number.  M  =  i  ^aj./cg,  is  assumed  to  be  small.  Axial  velocity  components  on 
the  left  and  right  boundaries  of  the  flow  configuration  are  assumed  to  be  strictly  equal 
to  the  shear  flow  velocity  In  other  words,  x'  =  0  and  x'  =  L'  represent  fixed 

nodal  surfaces  for  the  axial  acoustic  velocity.  This  configuration  is  highly  idealized, 
but  provides  an  opportunity  to  develop  an  analytical  investigation  of  trapped  acoustic 
signals  in  a  shear  flow,  when  multiple  wave  reflections  occur  on  the  axial  boundaries. 

The  dimensionless  equations  for  the  compressible,  viscous,  and  heat-conducting 
fluid  motion  in  the  duct  are  similar  to  those  in  Ref.  [5], 


p  = 

PT, 

(1) 
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=  0, 
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where  the  thermophysical  properties  have  been  assumed  constant  for  convenience. 
Subscripts  t,  x,  and  y  denote  partial  derivatives.  is  the  mean  flow 

Reynolds  number,  h  =  d' IV  the  aspect  ratio  of  the  duct,  Pr  the  Prandtl  number, 
and  the  nondimensional  dissipation  function.  Other  nondimensional  variables  are 
defined  in  terms  of  dimensional  quantities  by 


i  p.  P.  T) 


X  = 


£ 

V' 


iy'.p'.T) 
(Po-Po*  ^o) 
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V  = 
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(6) 


4 


The  subsequent  analysis  is  performed  in  the  parameter  ranges  relevant  to  how 
conditions  in  solid  rocket  engine  chambers  (excluding  the  exit  nozzle),  where  the 
values  of  M  and  R.  are  of  ^irid  (9(10*’).  respectively  61.  .\ccordingly.  it  is 

reasonable  to  define  the  perturbation  limit  as  1/ R.  —*  0;  M  —  0:  and  1//?,  C  M .  In 
addition,  one  assumes  that  the  aspect  ratio  h  <  0(1). 

Given  the  steady,  fully-developed  shear  flow  velocities  u  =  L'fy)  and  v  =  0.  one 
ob.serves  from  (l)-(5)  that  the  basic  pressure  distribution  is  independent  of  the  y 
direction,  and  that  dpidi  =  OiM^ / hRg),  negligible  compared  to  the  Oi\P)  acoustic 
pressure  to  be  considered.  The  high  Reynolds  number  limit  implies  that  viscous  and 
thermal  diffusion  exerts  little  influence  on  the  fluid  motion,  except  in  the  thin  acoustic 
boundary  layers  adjacent  to  solid  surfaces,  typically  of  thickness  less  than  1%  of  the 
duct  width  [5].  In  this  respect,  the  isentropic  relations, 

together  with  the  transport-free  version  of  (2)-{4)  (the  Euler  equations)  adequately 
describe  the  fundamental  physical  phenomena  in  the  core  flow  region. 

The  traveling  wave  study  in  Ref.  [5]  is  based  on  acoustic  velocity  disturbances 
of  0{M)  relative  to  the  shear  flow  velocity  U{y).  The  corresponding  acoustic  ther¬ 
modynamic  disturbances  must  be  0{\P).  a  magnitude  frequently  encountered  in 
stable  solid  rocket  motors,  where  an  absolute  pressure  oscillation  of  1  to  2%  is  often 
observed.  When  acoustic  perturbations  of  the  type 

u  =  U(y)  +  Mu,  u  =  Mv,  (S) 

p  =  1  +  ,V/^7p,  p  =  I -f  A/^p,  T  =  \  +  M^t,  (9) 

are  used  in  the  transport-free  version  of  (2)-(4),  the  following  equations  are  obtained 
for  the  acoustic  variables; 

p,  4-  4-  i-y  4-  .V/C'(y)pr  4-  [(pu)r  +  (pi))y]  =  0,  ( 10) 
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(11) 


lit  +  M  [U{y)u^  -t-  vL''{y)]  -r  .\r{uur  +  t'tly)  =  - 


Pr 


1  T  .\Pp' 


Vt  +  MC{y)Vr  +  M^iuVr  +  t-L’j,)  =  “777: 


Py 


(12) 


/i'll +.VP^]' 

Eqns.  (10)-(12)  show  that  the  mean  flow  i’(y)  affects  only  the  0{M)  acoustic  phe¬ 
nomena  in  the  duct.  The  quadratic  convective  terms,  representing  the  nonlinear  effect 
in  the  wave  system,  are  even  smaller,  of  0{\P). 

The  wave  field  evolves  from  an  initial  axially  distributed  velocity  disturbance: 

t  =  0,  u  =  Asmikirx),  v  =  p  =  Q.  (13) 

Since  the  left  and  right  boundaries  are  fixed  nodal  surfaces  for  u, 


X  =  0,  T,  u  =  0.  (14) 

The  normal  velocity  component  vanishes  on  the  impermeable  duct  sidewalls.  For 
mean  flow  Uiy)  symmetric  with  respect  to  the  center-plane  (y  =  0),  the  problem 
is  symmetric  and  can  be  solved  in  the  region  0  <  y  <  1  only.  The  appropriate  y 
boundaxy  conditions  are 

y  =  0, 1;  u  =  0.  (15) 


3  Solution  Development 

The  solution  to  the  lowest  order  approximation  to  the  equations  in  (10)-(15)  can  be 
written  as 


P\ 

=  —  ,Asin(fcxt)cos(fcTri), 

(16) 

=  Acos(ixf)sin(fc7ri), 

(IT) 

I'l 

=  0, 

(13) 

where  pi.  ui  and  tq  are  defined  by  the  asymptotic  expansions 

(p  =  ilTi -f.  A/'52  +  0(.V/^),  vp  =  (u,u,p,p.T).  (19) 
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It  follows  from  the  isentropic  relations  ''T).  i'3)  and  i'I9)  that 


Pi  =  Pi-  P2  =  P2.  (20) 

The  leading  order  solution  in  (16)-l  IS)  describes  an  axial  mode  with  a  harmon¬ 
ically  time-varying  amplitude,  independent  of  the  mean  flow.  In  order  to  obtain 
information  about  how  the  shear  flow  affects  the  acoustic  modes  through  convection 
and  refraction  effects,  the  next  order  solutions  must  be  constructed. 

The  0{M)  acoustic  equations  derived  from  (10)-(12)  and  (19)  take  the  form 


P2t  +  ^2r  A  t-'2y  —  (plPlri 

(21) 

U2t  +  P2i  =  -f’(y)uix  - 

(22) 

+  ^P2y  =  -f'(y)Vlx, 

(23) 

where  density  has  been  replaced  by  pressure  according  to  (20).  Eqns.  (16)-(18)  can 
be  combined  with  (21)-(23)  to  find  the  nonhomogeneous  wave  equation  for  p2, 

P7tt  -  +  ^P^yy)  =  -'2Ak'^rr^U{y)cos{kTrt)sm{k:rx).  (24) 

The  forcing  function  on  the  right  side  is  y-dependent  and  excites  fully  two-dimensional 
0[M)  acoustic  motion.  The  initial  conditions  for  pj  are  simply  (cf.  (13),  (19)  and 
(21)) 

p^[t  =  0)  =  P2tyt  =  0)  =  0.  (25) 

The  appropriate  boundary  conditions  are  derived  by  applying  the  results  in  (16)-(23) 
to  (U)-(15),  which  gives 

P2,(i  =  0)  =  —Ak-xC{y)cos{kwt)^ 

P2r(x=l)  =  -(-T)‘'.-4ArT6’(y)  cos(^Tf),  (26) 

P2y(y  =  0)  =  P:i,(y  =  1)  =  0.  (27) 


and 


Eqns.  (24)-(27)  constitute  a  well-defined  elementary  two-dimensional  hyperbolic  sys¬ 
tem  with  a  distributed  source  and  axial  boundary  excitations.  One  notices  that 
boundary  conditions  (26)  are  not  compatible  with  (25)  at  t  =  0.  .A.s  a  result,  a  pair 
of  propagating  discontinuities  in  pix,  or  weak  discontinuities  as  described  by  Landau 
and  Lifshitz  [7],  are  introduced  at  x  =  0  and  1  at  <  =  0. 

The  solution  to  (24)-(27)  can  be  written  aa 


—  =  —  cos(kTrt)  sinfArtrx)  '^Un  cosjmry) 

n=0  2 


+  ^  Y^emenC\rmnik,t)cos{rmrx)cos(niry),  (28) 

m=0(m^Ac)  n=0 

where  £m  and  £n  are  the  Neumann’s  number  whose  value  equals  1  if  the  subscript  is 


zero  and  2  if  it  is  a  positive  integer.  The  coefficients  and  r,n„ 


are  defined  by 


=2/;c'(,) 


cos{niry)dy, 


1  it  (m^  4- 

27r  (m^  +  fi^  -  (m^  — 


2k^  —  ,,  . 
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W?  ->rk^  ' 
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where  h  =  n/h.  If  m}  =  P,  T^n  should  be  evaluated  by  taking  the  limit 

4-  —  A:*  -♦  0  in  (30),  which  yields  0(f)  amplification  in  pj.  It  can  be  verified 

that  (28)  satisfies  equation  (24)  «is  well  as  all  the  initial  and  boundary  conditions. 

.A  more  concise  form  of  the  solution  is  derived  by  combining  the  two  terms  in  (28). 
which  yields 


where 

0rnn  — 


Ei 
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OO  PC  ^ 

H  £menOnQmn{k,t)cos{miTx)cos{niry), 
m=0(m#ik)  nsO 


(31) 


Jcos(A:7rf ) 


f  sin(fc7rf), 


4-  =  k^. 

(32) 
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In  the  derivation  of  (31),  sin(krx]  in  i2S)  is  expressed  in  terms  of  its  Fourier  cosine 
series  whose  derivative  converges  to  zero  instead  of  its  true  values  at  j  =  0  and  1. 
As  a  result  the  boundary  conditions  (26)  are  not  satisfied  by  (31).  However.  (31) 
gives  correct  x-derivatives  elsewhere  including  the  immediate  neighborhood  of  x  =  0 
and  1.  The  values  of  are  not  affected  at  all  because  the  series  given  by  (31)  is 
uniformly  convergent  throughout  the  entire  domain.  Thus  one  can  justifiably  use 
(  31)  to  evaluate  pi  everywhere,  and  its  first  derivatives  everywhere  except  at  x  =  0. 
1  and  at  the  locations  of  weak  discontinuities  where  these  derivatives  are  not  defined. 

The  Fourier  coefficients  of  the  double  series  in  (31)  are  proportional  to  Un  which 
represents  the  sheiir  flow  effects,  and  vary  with  k  and  t.  One  observes  from  (32)  that 
for  each  nonresonant  Fourier  mode,  the  coefficient  contains  two  harmonic  functions  of 
time  with  different  frequencies.  The  first  function,  cos{kirt),  has  a  frequency  equal  to 
that  of  the  distributed  forcing  function  in  (24),  or  that  of  the  imposed  leading  order 
acoustic  pressure.  It  describes  the  time  dependence  of  the  Helmholtz  mode  of  bulk 
oscillation  with  x  and  y-dependent  amplitudes. 

The  second  cosine  function  in  (32)  is  associated  with  wave  phenomena.  For  each 
pair  of  m  and  n,  the  product 

Pmn  —  cos  TTtJ  cos(m7ri)  cos(n7ry)  (33) 

represents  a  two-dimensional  oblique  mode,  which  can  be  decomposed  into  four  plane 
traveling  waves: 

P-nn  =  i  {cos  [>/m^  +  n‘iir[t-  cos  [v/m*  -h  ^  -  2+;)] 

-l-cos  [n/tti^  +  n^ir(^t-  x'jj)]  -I-  cos  -(-  tt  (<  -  x")] }  ,  (.34) 

where 

mx  ny  mi  -  ny 
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(35) 


—  mx  -t-  ny 
\/m^  -r 


—  mx  —  ny 


are  the  paths  taken  by  the  traveling  waves  reflecting  from  duct  walls  and  the  axial 
acoustic  reflectors.  The  wave  paths  are  oblique  in  general.  However,  when  either 
m  =  0  or  n  =  0,  the  traveling  waves  are  either  transverse  or  axial,  respectively, 
and  the  number  of  traveling  waves  reduces  to  two.  Since  (31)  contains  an  infinite 
number  of  Fourier  modes,  one  concludes  that  p-z  consists  of  infinite  numbers  of  axial, 
transverse  and  oblique  standing  waves,  in  addition  to  the  bulk  oscillation  mentioned 
earlier.  This  is  in  sharp  contr<ist  to  the  results  of  the  traveling  wave  study  [5]  where 
only  a  limited  number  of  wave  modes  can  propagate  along  a  long  duct.  The  multiple 
reflections  of  acoustic  signals  on  the  inserted  acoustic  reflectors  create  a  more  complex 
acoustic-shear  flow  interaction  response. 

Although  the  asymptotic  solution  given  by  (16),  (19)  and  (31)-(32)  is  derived 
primarily  for  processes  occurring  on  the  acoustic  time  scale  of  the  duct,  it  is  seen  to 
be  uniformly  valid  on  the  mean  flow  passage  time  scale  t  ~  0{M~^ )  as  well,  except  for 
the  resonant  C2ises.  Resonance  occurs  when  the  frequency  of  a  transverse  or  oblique 
standing  acoustic  wave  equads  to  that  of  the  forced  vibration.  It  causes  0(f)  amplitude 
growth  (cf.  (32))  and  invalidates  the  asymptotic  expansion  (19)  as  f  ~  0(A/“M- 
Comparing  to  the  0(f‘''*)  resonant  growth  found  for  refraction  induced  traveling 
waves  [5],  it  is  concluded  that  a  trapped  wave  system  produces  faster  growth  rate  for 
the  resonant  mode,  owing  to  the  effect  of  multiple  axial  reflections.  In  general,  the 
asymptotic  solution  remains  valid  until  f  ~  0(,V/~^),  the  time  required  for  the  small 
nonlinear  effects  to  accumulate,  which  leads  to  significant  deformation  of  the  leading 
order  waveform  and  possible  formation  of  a  weak  shock  in  the  confined  geometry  [8]. 

.\  complete  solution  should  include  the  thin  acoustic  bounadry  layers  adjacent 
to  the  duct  walls.  .An  approach  identical  to  that  used  in  Ref.  [5)  can  be  applied  to 
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reveal  the  structure  of  the  trausport-dominated  layers.  .A.  preliminary  analysis  shows 
a  basically  uniform  acoustic  pressure  distribution  across  the  laver,  whose  thickness  is 
proportional  to  the  ratio  of  the  fluid  kinematic  viscosity  to  the  mode  number  i  k]  of  the 
basic  core  pressure  disturbance.  .A  multitude  of  velocity  standing  modes  induced  by 
the  core  acoustic  flow  are  strongly  damped  by  viscous  and  thermal  diffusion  to  satisfy 
no-slip  and  appropriate  thermal  boundary  conditions  [5,  9].  Since  no  fundamentally 
new  discoveries  are  expected  from  a  detailed  analysis,  it  is  not  pursued  here. 

4  Example  Calculations  and  Discussion 

4.1  Acoustic  convection  and  refraction 

The  effect  of  a  sheared  mean  flow  on  acoustic  pro'-  o.  can  be  more  clearly  elucidated 
by  rewriting  (31)  as  p2  =  P2C  +  P7R  waere 

P2C  =  CmQm^jik*:roJ,mTrx),  (36) 

30  30 

P2fi  =  2.4  ^  V  c„C'n0mn  (fc.  Ocos{m7ri)  cos(ntry).  (37) 

i7i=0(m^*)  n=l 

The  first  component  of  p2,  p^c,  contains  all  the  n  =  0  Fourier  modes.  It  is  directly 
proportional  to  the  shear  flow  velocity  averaged  across  the  duct  {Uol2)  (cf.  (29)),  and 
is  the  correction  to  the  basic  acoustic  pressure  (16)  that  would  arise  in  a  pure  slug 
flow.  This  effect,  henceforth  referred  to  cis  the  bulk  convection  effect,  is  accountable 
for  the  axiad  standing  wave  disturbances  and  the  y-independent  part  of  the  bulk 
oscillations  in  p2.  The  magnitude  of  bulk  convection  increaises  with  the  average  shear 
flow  velocity  and  the  mode  number  (k)  of  the  basic  acoustic  disturbaince,  as  can  be 
shown  from  (32).  When  \fp2c  is  superimposed  on  the  baisic  acoustic  pressure  pt. 
the  total  acoustic  pressure  remains  one-dimensional,  although  its  waveform  is  slightly 
distorted  in  the  axial  direction  by  a  time-dependent  periodic  function. 
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The  quaxitity  p2fl  in  (37)  arises  from  convection  of  the  basic  axial  acoustic  dis 
turbance  at  the  local  differential  velocity  i'{y)  —  i’o/--  Thus,  is  associated  with 
refraction  of  the  axially  distributed  solution  pi  in  ( 16)  by  the  shear  flow.  It  is  seen 
to  be  intricately  related  to  the  specific  shear  flow  velocity  profile  through  the  Fourier 
coefficients  Un.  Both  transverse  and  oblique  standing  acoustic  waves  are  present,  and 
amplitude  amplification  of  certain  modes  is  possible  <is  a  result  of  resonance  phenom¬ 
ena.  The  P2R  solution  also  includes  bulk  oscillations  with  p-dependent  amplitudes  as 
mentioned  in  Section  3. 

The  general  features  of  the  acoustic  field  described  by  (16)  and  (31)-(32)  can  be 
revealed  graphically  by  evaluating  these  expressions  numericidly.  In  the  subsequent 
calculations,  each  Fourier  summation  is  carried  out  up  to  the  50th  term,  with  a 
truncation  error  of  less  than  10”^  based  on  comparisons  of  the  results  with  those 
from  summations  of  100  or  more  terms.  Three  types  of  shear  flows  f/(y),  as  well  as 
different  combinations  of  mode  number  k  and  duct  aspect  ratio  h,  are  employed  to 
generate  representative  results.  The  amplitude  of  the  initial  disturbance  A  is  kept  at 
unity  since  both  pi  and  pj  exhibit  the  same  simple  proportionality  to  it.  Results  are 
presented  in  terms  of  p2c  and  p^n,  mentioned  above,  for  easier  physical  interpretation. 

Figs  2a-2c  display  the  time-variations  of  the  three  acoustic  pressure  components 
on  the  duct  wall  (p  =  1),  at  axiaJ  locations  x  =  0,  1/4,  and  1/2,  respectively.  The 
wave  field  is  generated  by  the  relaxation  of  the  initial  velocity  disturbance  ti  =  sin(Tx) 
(^'  =  1)  in  a  fully  developed  laminar  flow  field  described  by  17  =  1  —  p^-  The  duct 
section  considered  has  equal  length  and  width,  so  that  h  =  1/2.  In  each  figure, 
the  solid  line  denotes  the  basic  acoustic  pressure  pi  as  calculated  from  (16).  The 
amplitude  of  its  oscillation  is  seen  to  decre«ise  from  maximum  at  the  antinode  x  =  0 
(cf.  Fig.  2a)  to  zero  at  the  nod<il  plane  x  =  1/2  (cf.  Fig.  2c).  The  long  and  short  dashed 
lines  represent  the  bulk  convective  correction  p^c  and  the  refractive  correction  pj/?, 
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calculated  from  (36)  and  (3i)  respectively.  The  former  is  a  (/-independent,  periodic 
function  of  time  whose  amplitude  and  structure  vary  with  logitudinal  position  x, 
as  observed  from  the  figures.  The  time  response  of  depends  on  both  spatial 
coordinates  and  exhibits  non-periodic,  irregular  oscillatory  behavior,  as  a  result  of 
the  superposition  of  signals  with  noncommensurate  frequencies  arising  from  transverse 
and  oblique  standing  modes  as  well  as  the  bulk  modes.  The  slope  discontinuities  on 
both  P2C  and  p^R  curves  denote  the  passage  of  the  weak  discontinuities  created  by 
the  initial  jump  in  p2x  (cf.  (26)).  The  curves  in  Figs.  2a-2c  are  plotted  up  to  t  =  10. 
at  which  the  mean  flow  completes  one  passage  through  the  duct  section  if  M  =  0.1. 
and  an  axial  acoustic  signal  is  reflected  10  times  on  the  left  and  right  boundaries. 
Longer-time  calculations  produce  similar  irregular  patterns  for  p2R  while  pi  and  p2c 
continue  their  periodic  variations.  It  is  particularly  interesting  to  notice  in  Fig.  2c 
that,  although  the  basic  acoustic  pressure  pi  vanishes,  the  total  acoustic  pressure 
p  =  Pi  -f  .V/(p2c  +  P2r)  is  nonzero  because  a  mean  flow  field  with  shear  is  present. 

The  two-dimensional  structure  of  the  acoustic  pressure  due  to  refraction  (p2ft) 
is  illustrated  in  Figs.  3a  and  3b,  which  depict  p2R  distributions  in  the  duct  for  an 
acoustic-shear  flow  system  characterized  by  A:  =  3,  A  =  1/6,  and  =  1  —  at 
t  =  5/8  and  7,  respectively.  The  two  snapshots  are  more  th«ui  6  acoustic  time  units 
apart,  and  thus  do  not  reflect  the  waveform  evolution  process,  which  occurs  on  a  much 
shorter  time  scale.  The  lines  of  sharp  slope  changes  cissociated  with  the  propagating 
weak  discontinuities  are  obvious  on  the  p2R  surfaces.  The  absolute  magnitude  of  the 
refraction  effect  is  smaJl  because  p2R  must  be  multiplied  by  M. 

4.2  Effect  of  mean  flow  type 

Figs.  4a  and  4b  illustrate  the  time  history  of  pi,  p2c  and  p2R  at  j  =  0,  y  =  1  under 
the  same  conditions  a.s  in  Fig.  2a,  except  for  different  shear  flow  velocity  types.  In 
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Fig.  4a,  U  =  1  —  y.  corresponding  to  a  fully-developed  Couette  flow  in  the  half 
duct  considered.  The  profile  U  —  (1  —  is  used  to  find  the  results  shown  in 
Fig.  4b.  It  is  easily  observed  from  Figs.  2a,  4a  and  4b  that  the  shear  flow  described 
by  the  one-seventh  power  law  generates  the  largest  bulk  convection  effect  p2c  because 
it  hats  the  largest  average  velocity.  Analogously  the  parabolic  flow  creates  a  larger 
bulk  convection  effect  than  the  linear  flow.  The  refraction  effects  for  the  linear  and 
parabolic  flow  types  are  of  comparable  magnitude.  The  one-seventh  power  law  flow 
generates  the  smallest  acoustic  refraction  because  the  velocity  gradient  is  relatively 
small  in  most  of  the  core  flow  region.  The  large  velocity  gradient  is  concentrated 
within  a  thin  layer  near  the  wail  which  is  too  narrow  to  promote  acoustic  refraction 
on  a  global  scale  in  the  core. 

The  refraction  induced  acoustic  pressure  fluctuations  at  specified  locations  are 
plotted  alone  in  Figs.  5a  and  5b  for  the  parabolic  and  one  seventh  power  mean  flow 
types,  under  conditions  it  =  1  and  /i  =  1/2.  The  solid  amd  dashed  lines,  representing 
P2r  at  an  axial  position  x  =  1/2  on  the  center-plane  {y  =  0)  and  duct  wall  (y  =  1) 
respectively,  resemble  each  other  but  are  completely  out  of  phaise,  implying  the  relative 
importance  of  refraction  in  the  duct.  It  is  strikingly  noteworthy  that  superpositions  of 
various  linear  wave  structures  resulting  from  shear  flow  interactions  with  a  simple  axial 
distanbance  can  produce  f2urly  irregular  pressure  signals,  which  might  be  mistakenly 
attributed  to  nonlinear  phenomena  had  the  data  been  collected  from  numerical  or 
experimental  investigations.  The  shear  velocity  described  by  the  one-seventh  power 
law  is  again  shown  to  generate  smaller  p^p  than  the  other  two  C2ises. 

4.3  Effect  of  mode  number  and  duct  aspect  ratio 

The  earlier  traveling  wave  study  [5]  suggests  that  the  magnitude  of  acoustic  refraction 
phenomena  is  controlled  by,  among  other  parameters,  the  ratio  of  the  wavelength  to 
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the  duct  width.  For  a  duct  of  fixed  width,  higher  frequency  waves  induce  larger 
refraction  effect.  In  the  present  analysis,  due  to  the  finite  axial  length  of  the  flow 
geometry,  the  relationship  between  poR  and  the  parameters  k  and  h  is  more  complex 
(cf.  (37)  and  (321)  .N'onetheless.  numerical  evaluations  show  that  when  k  is  away  from 
resonant  frequencies  of  the  duct  the  same  trend  is  apparent. 

The  nonresonant  p-2R  curves  in  Figs  6a  and  6b  are  plotted  at  j  =  0.  y  =  0 
and  1  using  the  same  shear  flow  velocity  U  =  1  —  The  k  and  h  values  are 
selected  such  that  while  k  increases  from  1  in  Fig.  6a  to  3  in  Fig.  6b.  the  aspect 
ratio  h  decretises  in  proportion  so  that  the  wavelength  to  duct-width  ratio  remains 
constant.  .As  expected,  the  refraction  induced  pressure  fluctuations  in  the  two  figures 
are  comparable  in  magnitude,  and  the  frequency  of  fluctuation  increases  with  k.  The 
p2R  signals  are  all  fairly  irregular.  The  location  x  =  0  is  chosen  in  Figs.  6a  and  6b 
because  it  corresponds  to  the  anti-node  for  pi,  where  the  p^R  curves  are  of  typical 
magnitude  for  both  cases.  At  other  i-positions  p-zR  fluctuates  in  a  similar  manner 
with  different  rates  and  amplitudes  (see,  for  example.  Fig.  5a). 

Fig.  7  shows  a  near  resonant  case  that  results  in  a  long  period  beat  pattern.  Here 
A:  =  4,  =  1/4,  so  that  for  m  =  1  and  n  =  1,  -I-  fi*  =  17  while  k^  —  16. 

The  pzR  curves  are  plotted  at  x  =  1/8,  y  =  0  and  1,  for  (7  =  1  —  y^.  Beats  of 
even  longer  periods  and  leirger  maximum  amplitudes  can  occur  at  suitable  parameter 
combinations.  If  the  numerical  evaluation  of  the  solution  had  been  limited  to  times 
less  than  half  the  beat  period,  the  result  could  have  been  erroneously  interpreted  as 
an  amplified  mode. 

Fully  resonant  solution  behavior  is  shown  in  Figs.  8a  and  8b,  where  pzR  on  the 
duct  wall  is  depicted  at  x  =  1/10  and  x  =  1/2,  respectively,  for  =  1  —  y^.  In  Fig.  Sa 
k  =  0,  h  =  1/5.  so  that  the  transverse  mode  (m  =  0,  n  =  1)  is  amplified.  Fig.  8b 
corresponds  to  the  condition  k  =  o  and  h  =  1/ \/^.  which  promotes  the  amplification 


of  an  oblique  mode  (m  =  2,  n  =  1).  One  may  observe  beat  patterns  superimposed 
on  the  long-time  growth  arising  from  the  interference  by  near-resonant  modes. 

The  considerable  enhancement  of  the  refractive  acoustic  pressure  by  resonance 
requires  a  new  asymptotic  solution  on  the  flow  time  scale  t  ~  0{M~^ )  to  accomondate 
refraction  phenomena  in  the  leading  order  acoustics.  Likewise,  when  certain  modes 
oscillate  at  near-resonant  frequencies  of  the  system,  or  when  the  mode  number  k 
becomes  large  and  the  aspect  ratio  h  small,  so  that  .V/pj  ~  0(pi),  a  new  theory  needs 
to  be  developed. 

4.4  Inferences  about  solid  rocket  engine  acoustics 

Solid  rocket  combustion  chamber  dimensions  are  usually  associated  with  radius  to 
length  ratios  less  than  or  equal  to  0(1).  Experimental  observations  suggest  that  the 
first  few  acoustic  modes,  especially  the  k  —  I  mode,  are  most  frequently  encountered 
in  rocket  chambers.  Based  on  the  above  analysis,  the  convection  and  refraction  in¬ 
duced  acoustic  pressure  correction  terms  are  at  most  O(.Vf)  relative  to  the  leading 
order  acoustic  pressure,  or  0{M^)  relative  to  the  overadl  pressure.  Although  reso¬ 
nance  in  P2fi  significantly  magnifies  the  acoustic  refraction,  it  occurs  rarely  for  low 
k  modes.  In  fact  the  resonant  condition  rn}  -f  (nfhy  =  k^  indicates  that  there  can 
be  no  resonant  modes  \{  kh  <  1,  which  is  likely  to  be  the  case  given  the  sizes  of  k 
and  h  mentioned  above.  Thus,  a  mean  flow  with  .V/  =  0.1  causes  at  most  O(l0~M 
changes  in  acoustic  pressure,  or  0(10“^)  changes  in  total  pressure.  The  nonlinear 
effects,  represented  by  quadratic  terms  in  the  convective  operators  and  responsible 
for  acoustic  streaming,  are  0(10”*)  smaller.  They  are  of  no  importance  on  the  time 
scale  considered  here  and  accumulate  only  on  the  longer  time  scaJe  t  ~  O(IO^). 

The  resonant  phenomenon  discussed  above  appears  in  second  order  acoustic  quan¬ 
tities  only,  and  is  excited  by  interactions  with  a  sheared  mean  flow.  It  should  be 
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distinguished  from  the  more  powerful  resonance  of  leading  order  acoustic  pressure  in 
a  rocket  engine  configurations,  often  associated  with  excitations  arising  from  chemical 
heat  release,  which  is  beyond  the  scope  of  this  paper. 

5  Summary 

In  this  paper  a  mathematical  model  is  established  to  study  the  interaction  between 
a  low  Mach  numbei  parallel  shear  flow  and  an  a.xially  distributed  acoustic  distur¬ 
bance  trapped  in  a  section  of  a  planar  duct.  The  solution  development  is  based  on 
asymptotic  expansions  defined  for  the  small  mean  flow  Mach  number  limit  {M  -<■  0). 
The  Fourier  series-based  analytical  solution  in  terms  of  standing  modes  gives  explicit 
modal  response  of  the  acoustic-shear  flow  interactions.  The  major  conclusions  are: 

•  Bulk  (slug  flow)  convection  induces  purely  axial  standing  acoustic  waves  and 
^-independent  bulk  oscillations  at  the  frequency  of  the  imposed  initial  distur¬ 
bance.  The  magnitude  of  the  bulk  convection  effect  increases  with  the  average 
shear  flow  velocity  and  the  mode  number  {k)  of  the  leading  order  acoustic  pres¬ 
sure. 

•  Refraction  generates  purely  transverse  and  oblique  standing  waves  as  well  as 
bulk  oscillations  with  both  x-  and  j/-dependent  amplitudes.  The  magnitude 
of  refraction  increases  with  the  shear  flow  velocity  gradient,  the  duct-width  to 
length  ratio  (A),  and  the  mode  number  of  the  leading  order  acoustic  pressure 
ik).  The  refraction  size  also  increases  if  the  frequencies  of  the  induced  acoustic 
modes  become  close  to  the  resonant  frequency. 

•  Combination  of  the  above  effects  gives  irregular  pressure  signals  at  given  po¬ 
sitions  in  the  duct.  They  are  of  0{.V/)  relative  to  the  imposed  leading  order 
acoustic  pressure. 
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•  Axial  wave  reflections  cause  faster  resonant  growth  of  acoustic  refraction.  Res¬ 
onance  does  not  occur  for  hk  <  1. 

•  Refraction  effects  are  very  small  in  both  traveling  and  standing  wave  model 
studies,  in  the  parameter  ranges  of  solid  rocket  engines. 

The  primary  point  is  that  quite  elementciry  disturbances  can  evolve  into  surpris¬ 
ingly  complex  wave  structures  as  a  result  of  refraction  effects,  although  these  effects 
are  small  in  the  parameter  ranges  of  solid  rocket  engines.  This  simplified  analysis 
gives  insight  into  the  acoustic  modal  interactions  with  the  shear  flow,  which  have  not 
been  extracted  from  computational  data  in  past  studies.  In  fact  there  appears  to  be 
a  specific  need  to  develop  data  analysis  tools  that  are  capable  of  using  grid  point  flow 
data  to  determine  the  propagating  wave  patterns  present  in  an  acoustically  excited 
shear  flow. 
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Figure 

1.  Schematic  of  a  parallel  shear  flow  through  a  planar  duct  section  containing 
acoustic  disturbances. 

2.  Time  variations  of  pi,  p^c  and  P2R,  as  calculated  from  (16),  (36)  and  (37) 

respectively,  on  the  duct  wadi  y  =  1  at  three  axial  locations;  (a)  i  =  0;  (b) 

X  =  1/4;  (c)  X  =  1/2.  The  acoustic  field  is  generated  by  the  relaxation  of  an 
initial  disturbance  u{t  =  0)  =  sin(Tx)  in  a  laminar  duct  flow  U  =  I  —  y^.  The 
duct  aspect  ratio  h  ■=  1/2. 

3.  Distribution  of  second  order  acoustic  pressure  due  to  refraction  (p^y?)  in  the 

duct  section  at  (a)  t  =  5/S;  (b)  t  =  7,  as  calculated  from  (37).  The  refraction 

effects  are  generated  as  the  initial  acoustic  disturbance  u(t  =  0)  =  sin(37rxj 
relaxes  in  the  shear  flow  field  U  =  1  —  The  duct  aspect  ratio  h  =  1/6. 

4.  Time  variations  of  p\,  p-ic  and  p-iR  at  x  =  0,  y  =  1,  generated  by  the  relaxation 
of  the  initial  disturbance  u{t  =  0)  =  sin(Tx)  in  shear  flow  fields  described  by 
(a)  T'  =  1  —  y;  (b)  u  =  (1  —  y)'*^^.  The  duct  eupect  ratio  h  =  1/2. 

5.  Time  variations  of  refractive  acoustic  pressure  p2R  as  calculated  from  (37),  at 
locations  x  =  1/2,  y  =  0  and  1.  The  refraction  effects  are  generated  by  the 
relaxation  of  the  initial  acoustic  disturbance  u{t  =  0)  =  sin(ffi)  in  shear  flows 
(a)  U  =  I  —  y^;  (h)  U  =  (1  —  y)*^^.  The  duct  aspect  ratio  h  =  1/2. 

6.  Time  variations  of  refractive  acoustic  pressure  p2R  at  locations  x  =  0,  y  =  0 
and  1,  generated  as  the  initial  acoustic  disturbance  u(t  =  0)  =  sin(^'Tx)  relaxes 
in  the  shear  flow  field  U  =  I  —  y^.  (a.)  k  =  1.  h  =  1/2;  (b)  k  —  3.  h  =  1/6. 
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7.  Near- resonant  variations  of  refractive  acoustic  pressure  p2n  with  time  at  loca¬ 
tions  X  =  1,  3.  7  =  0  and  1.  as  the  initial  acoustic  disturbance  ui.t  =  0)  = 
sin(4Tx)  relaxes  in  the  shear  flow  field  6'  =  1  -  i/.  The  duct  aspect  ratio 
h  =  1/4. 

3.  Resonant  amplifications  of  refractive  acoustic  pressure  p2R  on  the  duct  wall 
(y  =  1),  at  specified  axial  locations,  as  the  initial  acoustic  disturbance  iiit  = 
0)  =  sin(^7rx)  relaxes  in  the  shear  flow  field  U  =  I  —  y^.  (a)  k  =  o,  h  =  1/5. 
X  =  1/TO:  (b)  k  =  0,  h  =  x  =  1/2. 
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